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It  is  a  well  known  fact  that  high  frequency,  or  radio  frequency,  radiation  can  be  directly  harmful  to  biological 
tissue.  The  radiation  frequencies  to  which  humans  are  most  exposed,  however,  are  the  extremely  low  frequencies, 
or  ELFs.  To  date  there  has  been  no  definitive  measure  for  predicting  the  ability  of  low  frequency  electromagnetic 
radiation  to  cause  adverse  biological  efiects.  Two  specific  measures,  average  magnetic  field  and  kinetic  index, 
have  been  used  in  studies  to  determine  if  they  are  usefiil  predictors  of  adverse  effects,  specifically  leukemia. 
Unfortunately,  both  have  shown  only  marginal  results. 

This  thesis  analyzed  and  improved  the  kinetic  index  measure  (Thomas,  et  al.,  1994)  by  utilizing  laboratory' 
biological  data  to  increase  the  validity  of  the  parameters  of  an  existing  model.  Then,  using  this  “improved” 
model,  raw  magnetic  field  exposure  data  from  a  previous  study  was  analyzed  and  compared  with  the  measurement 
method  used  in  that  study  (average  magnetic  field). 

The  results  of  the  model  modification  showed  that  accurate  modeling  of  the  existing  laboratory  data  could 
only  be  accomplished  by  including  and  additional  parameter  to  the  existing  Litovitz  multistage  model.  The 
results  of  the  model  application,  however,  were  not  particularly  conclusive.  The  case-control  study  used  for  the 
model  qjplicatioo  was  less  than  ideal  requiring  that  more  rigorous  epidemiological  studies  be  conducted  in  order 
ID  aocarately  test  the  improved  model. 
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Abstract 


It  is  a  well  known  fact  that  high  frequency,  or  radio  frequency,  radiation  can  be 
directly  harmful  to  biological  tissue.  The  radiation  frequencies  to  which  humans  are  most 
exposed,  however,  are  the  extremely  low  frequencies,  or  ELFs.  To  date  there  has  been  no 
definitive  measure  for  predicting  the  ability  of  low  frequency  electromagnetic  radiation  to 
cause  adverse  biological  effects.  Two  specific  measures,  average  magnetic  field  and 
kinetic  index,  have  been  used  in  studies  to  determine  if  they  are  useful  predictors  of 
adverse  effects,  specifically  leukemia.  Unfortunately,  both  have  shown  only  marginal 
results. 

This  thesis  analyzed  and  improved  the  kinetic  index  measure  (Thomas,  et  al., 

1994)  by  utilizing  laboratory  biological  data  to  increase  the  validity  of  the  parameters  of 
an  existing  model.  Then,  using  this  “improved”  model,  raw  magnetic  field  exposure  data 
from  a  previous  study  was  analyzed  and  compared  with  the  measurement  method  used  in 
that  study  (average  magnetic  field). 

The  results  of  the  model  modification  showed  that  accurate  modeling  of  the 
existing  laboratory  data  could  only  be  accomplished  by  including  and  additional 
parameter  to  the  existing  Litovitz  multistage  model.  The  results  of  the  model  application, 
however,  were  not  particularly  conclusive.  The  case-control  study  used  for  the  model 
application  was  less  than  ideal,  requiring  that  more  rigorous  epidemiological  studies  be 
conducted  in  order  to  accurately  test  the  improved  model. 
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1.  Introduction 


1.1  Background 

There  is  little  doubt  within  the  scientific  community  that  extremely  low  frequency 
(ELF)  electromagnetic  fields  (EMF)  can  produce  some  types  of  biological  changes.  What 
is  in  dispute  is  what  types  of  effects  these  EMFs  have  and  whether  or  not  they  have  the 
potential  to  adversely  affect  human  health. 

One  type  of  model  which  has  been  proposed  and  attempts  to  address  magnetic 
field-induced  changes  at  the  cellular  level  is  the  Litovitz  multistage  DNA  transcription 
model  (Litovitz,  et  ah,  1990:297-312;  Litovitz,  et  ah,  1992:237-246).  This  model  seeks 
to  illustrate  the  effects  of  a  low-intensity,  pulsed,  ELF  on  the  transcription  of  messenger 
ribonucleic  acid  (mRNA)  in  terms  of  a  three-stage  process.  Its  prediction  asserts  that 
when  a  magnetic  field  is  turned  on  and  subsequently  held  at  a  constant  level,  the  mRNA 
concentration  initially  increases  and  then  decreases  to  some  basal  level.  The  rise-time  of 
the  peak  concentration  and  the  resulting  peak  concentration  of  mRNA  both  depend  upon 
the  initial  field  strength.  If  the  field  is  changed,  the  modeled  mRNA  concentration  can 
take  a  much  longer  amount  of  time  to  return  to  its  basal  level  than  if  the  field  had 
remained  energized  at  its  previous  level. 

Litovitz  ,et  ah,  chose  their  model  parameters  based  solely  upon  characteristics 
(not  specific  data)  exhibited  by  the  snapshot  responses  of  several  earlier  laboratory 
experiments  (Goodman,  et  ah,  1983:1283-1285;  Goodman,  et  ah,  1986:23-29).  These 
experiments  measured  timed  mRNA  responses  after  being  exposed  to  one  of  several 
different  types  of  electromagnetic  waveforms.  Litovitz  also  assumed  pulsed  EMFs  for  his 
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model  rather  than  EMFs  which  were  ubiquitous  for  the  general  population  (i.e.  60  Hz 
continuous).  Very  recent  research  (Lin,  et  al.,  1995),  however,  exposed  cells  to 
continuous,  interrupted,  and  changing  60  Hz  fields.  This  data  is  to  be  used  in  this 
research  to  improve  upon  the  Litovitz  model. 

1.2  Problem,  Objectives,  and  Scope 

The  major  problem  with  ELF-EMF  exposure  is  that  to  date  there  has  been  no  solid 
statistically  plausible  antecedent-consequence  relationship  drawn  between  these 
exposures  and  incidents  of  leukemia.  Upon  reading  Chapter  3  one  will  see  that  the  only 
developed  antecedents  that  have  maintained  a  statistical  correlation  with  ELF-EMF 
exposure  are  the  Wertheimer-Leeper  wiring  codes.  These  codes  assign  exposure 
categories  based  upon  the  wiring  configurations,  mainly  high-voltage  distribution, 
surrounding  the  individual  or  population  being  tested.  The  coding  is  a  partially  subjective 
process,  and,  therefore,  scientists  are  attempting  to  find  a  measure  which  will  both 
provide  good  correlation  and  be  objective. 

One  objective  measure  which  would  be  the  obvious  first-choice  is  the  average 
magnetic  field  exposure.  After  all,  this,  along  with  incidents  of  leukemia,  can  be 
measured  directly  and  would  intuitively  seem  to  be  the  best  choice.  London,  et  al.,  in 
1994  conducted  a  case-control  study  (Section  3.2  and  Chapter  5)  hoping  to  show  just  such 
an  antecedent-consequence  relationship.  Although  the  point  estimates  for  the  odds  ratios 
were  higher,  the  confidence  intervals  were  too  wide  to  give  this  the  required  level  of 
statistical  assurance. 
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The  Litovitz  multistage  model  is  the  basis  for  another  candidate  objective 
measure.  The  drawback  to  this  measure,  however,  is  two-fold.  First,  and  foremost,  the 
model  parameters  used  in  creating  this  measure  were  not  based  on  any  specific  biological 
data.  They  were  chosen  to  show  the  subjective  characteristics  which  had  been  discovered 
in  earlier  cellular  experiments.  In  the  one  study  which  was  conducted  using  this  measure, 
model  and  its  initial  parameters  as  the  basis  for  an  exposure  metric  (Thomas,  et  ah,  1994), 
the  researchers  were  not  able  to  show  a  correlation  between  the  calculated  metric  and 
incidents  of  childhood  leukemia  (section  3.6). 

This  research  will  propose  yet  another  possible  predictor  for  leukemia  incidents. 
This  measure  will  be  developed  by  modifying  the  existing  Litovitz  model  parameters 
including  the  possible  addition  of  parameters,  in  order  to  better  represent  mRNA 
responses.  This  will  be  accomplished  by  using  the  biological  data  from  Lin,  et  al.,  1995. 
Upon  developing  the  modified  model  we  will  then  use  this  model  to  predict  mRNA 
synthesis  from  a  data  set  taken  of  EMF-exposed  individuals  in  the  workplace.  Finally,  we 
will  use  these  predictions  to  calculate  real  world  occupational  risk  analysis  figures  which 
will  hopefully  show  an  improved  antecedent-consequence  relationship  between  ELF- 
EMF  and  incidents  of  leukemia  in  certain  occupational  groupings. 

The  potential  connection  that  this  model  has  with  cancer  formation  is  in  the 
altered  mRNA  production.  Cancer,  in  its  most  simplistic  definition,  is  the  uncontrolled 
reproduction  of  abnormal  cells.  As  the  next  chapter  will  point  out,  mRNA  plays  a  pivotal 
role  in  cell  reproduction.  If  the  Litovitz  model  is  correct,  and  mRNA  concentration 
rapidly  increases  upon  exposure  to  ELF-EMFs,  then  it  is  entirely  possible  that  these  fields 
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could,  at  the  very  least,  accelerate  the  reproduction  of  existing  abnormal  cells,  and  may 
actually  contribute  to  the  formation  of  abnormal  cells. 

1.3  Assumptions 

There  are  three  main  assumptions  that  deal  directly  with  the  Litovitz  multistage 
model  and  are  also  foundations  of  the  modified  model.  First,  “there  is  some  kind  of 
biological  transient  response  to  an  impressed  EMF-ELF  field.”  Second,  “there  is  a 
maximal  response  that  increases  proportionally  with  the  strength  of  the  irradiating  field.” 
Finally,  “there  is  a  peak  location  that  shifts  to  an  earlier  time  as  the  field  strength  is 
increased”  (Litovitz,  et  al.,  1990:298).  Additional  assumptions  proposed  by  Litovitz,  et 
al.,  and  maintained  throughout  the  course  of  this  research  are: 

1.  Exposure  to  an  EMF-ELF  will  increase  both  the  synthesis  and 
degradation  rates  in  the  multistep  model  (Figure  3  and  Figure  4)  (Litovitz, 
et  al.,  1992:241). 

2.  Reactions  (Figure  3  and  Figure  4)  are  strongly  biased  in  the  forward 
direction,  and  that  only  these  rate  constants  (with  the  addition  of  any 
parameters  added  by  this  research)  need  be  considered  (Litovitz,  et  al., 

1990:299). 

3.  The  nucleotide  reservoir  will  be  rapidly  replenished  ([A]  in  Figure  4) 

(Litovitz,  et  al.,  1990:300).  [This  means  that  the  supply  of  nucleotides 
with  which  to  maintain  mRNA  synthesis  will  be  limitless]. 
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2.  Fundamentals 


This  chapter  presents  a  review  of  the  principles  and  characteristics  of  magnetic 
fields  and  a  few  proposed  and  established  biological  effects  exhibited  by  them.  It  will 
also  describe  deoxyribonucleic  acid  (DNA)  and  mRNA  as  well  as  the  role  mRNA  plays 
in  transcribing  DNA  protein.  In  understanding  these  concepts  one  should  be  able  to 
follow  the  background  and  development  of  the  modified  model  which  this  thesis 
addresses. 

2.1  Magnetic  Fields 

While  electric  fields  are  created  by  the  existence  of  electric  charges,  magnetic 
fields  exist  only  when  those  charges  are  moving.  The  force  exerted  on  a  moving  charge 
by  a  magnetic  field  is  given  by  the  equation: 

F  =  q(vxB)  (1) 

where  F  is  the  force,  q  is  a  unit  charge,  v  is  the  magnitude  and  direction  of  the  relative 
motion  of  the  field  and  B  is  the  magnetic  flux  density.  The  unit  of  measure  which  will  be 
used  to  describe  the  magnetic  flux  density  is  the  milligauss  (mG). 

Normally  it  is  unusual  to  find  references  which  single  out  magnetic  fields.  They 
are  typically  referred  to  in  combination  with  electric  fields  in  a  quantity  known  as 
electromagnetic  radiation.  In  most  situations  these  two  fields  travel  in  tandem  and  are 
considered  as  a  single  entity  when  studying  the  effects  on  living  tissue.  In  the  higher 
frequency  ranges  such  as  those  involving  radio,  television,  light,  and  microwave  ovens 
this  is  the  case.  It  is  also  the  case  that  at  these  higher  frequencies  electromagnetic 
radiation  has  a  thermal  effect  on  biological  tissue.  However,  in  the  lower  frequency 
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ranges,  such  as  those  associated  with  the  transmission  and  distribution  of  electrical 
power,  electric  fields  and  magnetic  fields  can,  and  almost  always  are,  considered 
separately  and  do  not  outwardly  appear  to  have  an  effect  on  the  human  body.  The 
distinction  between  these  two  frequency  groups  is  well  established:  “When  the  distance 
from  a  source  of  electromagnetic  energy  is  large  compared  to  the  wavelength,  the  electric 
and  magnetic  fields  are  linked  and  are  considered  together  as  a  single  field.  When  the 
distance  from  the  source  is  small  with  respect  to  the  wavelength,  the  electric  and 
magnetic  fields  are  not  linked  and  can  be  considered  as  separate  entities”  (EPRI,  1989:4). 

The  separation  of  these  fields  at  low  frequencies  allows  the  investigation  of  the 
individual  effects  of  exposure  from  each  component  field.  The  consensus  of  the  scientific 
community  is  that  most  biological  effects  which  are  allegedly  associated  with  extremely 
low  frequency  (ELF)  electromagnetic  fields  (EMF)  are  most  likely  due  to  the  magnetic 
component  (Hileman,  1993:18).  The  reason  for  this  is  two-fold.  First,  trees,  building 
materials,  and  clothing  inadvertently  shield  the  human  body  from  most  electrical  field 
exposures.  Second,  the  body  itself  is  able  to  attenuate  electric  fields  significantly  due  to 
its  composition.  The  field  transmitted  to  the  inside  of  the  body  as  a  result  of  this 
attenuation  tends  only  to  be  one  millionth  to  one  hundred  millionth  of  the  applied  field 
(Hileman,  1993:18). 

Magnetic  fields,  on  the  other  hand,  pass  through  biological  material  unimpeded. 
This  occurs  due  to  the  lack  of  magnetic  materials  in  living  organisms.  Additionally,  any 
voltage  gradients  created  across  cell  walls  by  magnetic  fields  are  insignificant  compared 
to  those  created  by  the  cells’  endogenous  thermal  fields.  The  highest  voltage  gradient 
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which  can  be  initiated  by  a  60  Hz  magnetic  field  (1,000  to  5,000  mG)  across  a  cell  wall  is 
in  the  neighborhood  of  50  volts  per  meter.  The  natural  gradients  across  those  same  cell 
walls  averages  lO’  volts  per  meter  (Hileman,  1993:19).  Finally,  magnetic  fields  do  not 
have  enough  energy  to  heat  tissue  or  disrupt  DNA  directly  (Hileman,  1993:25). 
Considering  the  lack  of  magnetic  material  in  the  human  body,  the  relatively  insignificant 
voltage  gradients  created  by  an  exogenous  field,  and  the  low  energy  of  magnetic  fields,  it 
would  be  logical  to  conclude  that  magnetic  fields  should  have  little  or  no  effect  on 
biological  systems. 

Nevertheless,  they  have  been  shown  to  affect  the  synthesis  rate  of  DNA  through 
RNA  transcription  (Hileman,  1993:25).  This  phenomena,  which  will  be  described  in 
detail  later,  has  been  studied  and  modeled  by  Litovitz,  Montrose,  et.  al.  (Litovitz  et  al., 
1990:297-3 12)  and  forms  a  partial  basis  for  the  formulation  of  this  paper’s  proposed 
model. 

2.2  Natural  and  Typical  Exposures 

All  living  things  are  exposed  to  a  certain  amount  of  constant  electromagnetic 
energy  which  has  been  present  ever  since  the  earth  was  formed.  One  type  is  the  static 
field  which  exists  between  the  earth  and  the  ionosphere  (essentially  a  spherical  capacitor). 
The  most  visible  evidence  of  its  existence  is  the  occurrence  of  lightning  which  serves  to 
equalize  the  charges.  This  field,  however,  is  largely  electric. 

The  other  natural  field  which  occurs  is  the  geomagnetic  field  which  is  evident  in 
the  use  of  compasses.  Like  the  capacitive  field,  it  is  largely  static,  oscillating  slowly  over 
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a  24-hour  period,  but  in  contrast  it  is  largely  magnetic.  Its  strength  varies  from  300  to 
600  mG  between  the  equator  and  the  magnetic  poles  (Newman,  1992:1715). 

Man-made  sources  of  electromagnetic  energy  are  numerous.  High  voltage  power 
lines  which  receive  the  most  scrutiny  and  are  the  most  readily  visible  can  produce 
magnetic  fields  of  about  350  mG. 

In  the  home,  there  are  several  sources  of  magnetic  fields  ranging  from  hair  dryers 
to  electric  stoves.  Fields  in  close  proximity  to  these  objects  can  be  significantly  higher 
than  those  found  near  high  voltage  lines.  Table  1  below  shows  typical  values  of  magnetic 
fields  near  several  typical  household  appliances. 


Table 

Typical  Household  Magnetic  Exposures 

Appliance 

Magnetic  Flux  Density,  mG 

3  cm  30  cm  1  m 

10000-20000 

35-300 

0.7-10 

60-20000 

0.1-70 

<0.1-3 

Electric  shavers 

150-15000 

0.8-90 

<0.1-3 

Drills 

4000-8000 

20-35 

0.8-2 

Mixers 

600-7000 

6-100 

0.2-2.5 

Portable  Heaters 

100-1800 

1.5-50 

0. 1-2.5 

Blenders 

250-1300 

6-20 

0.3-1.2 

Television 

25-500 

0.4-20 

0.1-1.5 

Irons 

80-300 

1.2-3 

0.1-0.25 

Coffee  makers 

18-250 

0.8-1.5 

<0.1 

5-17 

0.1 -'2.5 

<0.1 

(EPRI,  1989:7) 


Electrical  distribution  lines  provide  three  specific  sources  of  magnetic  fields. 
Electrical  currents  in  the  primary  feeders  to  service  transformers  create  magnetic  fields  as 
do  the  secondary  service  drops  from  those  transformers  to  residential  or  commercial 
customers.  What  is  not  normally  considered,  however,  is  the  “net”  field  created  by  the 
vector  sum  of  the  secondary,  primary,  and  neutral  currents.  This  field’s  strength  is 
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inversely  proportional  to  the  distance  from  the  aggregate  source  while  the  field  strengths 
of  the  primaries  and  secondaries  are  inversely  proportionally  to  the  square  of  the  distance 
from  their  respective  sources.  This  net  field  is  not  only  stronger,  but  more  uniform  than 
those  from  its  component  fields. 

Finally,  many  homes  and  businesses  have  electrical  systems  which  are  grounded 
to  water  pipes.  This  practice,  though  safe,  produces  non-uniform  magnetic  fields  due  to 
the  return  neutral  current  from  home  appliances.  It  also  produces  fields  which  vary 
greatly  in  time  since  the  ground  current  changes  every  time  a  120-V  appliance  is  switched 
on  or  off  (EPRI,  1989:8).  In  several  studies  investigating  the  effects  of  magnetic  and 
electric  fields  on  humans,  the  Midwest  Research  Institute  (MRI)  in  Kansas  City,  Missouri 
consistently  found  that  the  greatest  effects  occurred  just  after  the  field  was  turned  off  or 
on  (Hileman,  1993:24). 

2.3  DNA,  mRNA,  and  Protein  Synthesis 

DNA  is  the  basic  genetic  code  of  which  all  living  things  are  constructed.  It 
consists  of  two  complementary  chemical  strands  which  are  coiled  around  each  other  to 
form  a  double  helix  structure.  The  sequenced  chemicals  (or  bases)  making  up  the  strands 
are  called  nucleotides.  In  DNA,  nucleotides  consist  of  a  sugar  and  a  phosphate  group 
plus  either  of  two  purine  bases  adenine  (A)  and  guanine  (G)  or  either  of  two  pyrimidine 
bases  thymine  (T)  and  cytosine  (C).  RNA  is  made  up  of  the  same  basic  units  except  that 
instead  of  T,  it  utilizes  the  pyrimidine  uracil  (U)  (CIE,  1994). 

The  two  sides  of  the  DNA’s  double  helix  structure  are  held  together  by  hydrogen 
bonds.  Each  base  on  a  strand  of  DNA  pairs  only  with  its  complement  on  the  other  strand. 
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G  pairs  only  with  C,  and  A  pairs  only  with  T  (or  with  U  in  RNA).  Each  three-base  set  on 
a  strand  codes  for  a  specific  amino  acid.  The  sequencing  of  these  amino  acids  is  the  basis 
for  protein  synthesis  (CEE,  1994). 

Protein  synthesis  occurs  in  ribosomes  which  are  located  in  the  cytoplasm.  DNA  is 
located  in  the  nucleus.  Data  needed  for  protein  construction  are  transferred  from  the 
nucleus  t  the  ribosomes  via  mRNA.  mRNA  is  closely  related  to  DNA  and  can  carry 
genetic  messages.  In  the  transcription  process  DNA  unwinds  and  separates  its  strands  so 
that  complementary  strands  of  mRNA  can  be  assembled  on  them.  A  strand  of  mRNA 
which  has  duplicated  the  DNA  information  then  travels  out  of  the  nucleus  to  the 
ribosomes,  where  transcription  of  the  actual  DNA  and  protein  synthesis  begin.  Since 
code  transmission  from  DNA  to  mRNA  is  extremely  Critical,  any  error  in  the  code 
ultimately  affects  protein  synthesis.  If  the  error  is  serious  enough,  it  eventually  affects 
some  body  trait  or  feature  of  the  subject  organism. 
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3.  Supporting  Research  and  Model  Development 


The  purpose  of  this  chapter  is  to  familiarize  the  reader  with  the  existing  biological 
research  into  the  EMF-leukemia  link,  the  EMF-biological  systems  link,  the  Litovitz 
multistage  model  which  characterizes  the  biological  link  mathematically,  and  a  study 
which  has  incorporated  this  model.  Through  the  epidemiological  studies  this  material 
should  present  a  clear  picture  of  the  concerns  EMFs  present  to  the  medical  and  scientific 
communities.  It  should  next  give  sufficient  background  to  allow  comprehension  of  the 
laboratory  experiments  which  have  sought  to  directly  link  EMFs  to  specific  biological 
effects.  And  finally,  one  should  be  able  to  understand  the  motivation  and  development  of 
the  current  Litovitz  multistage  toxicokinetic  model. 

3. 1  Residential  Epidemioiogical  Research 

The  most  cited  studies  which  cover  the  association  of  ELF-EMF  exposures  to 
cancer  are  a  group  of  three  childhood  case-control  studies,  two  in  Denver,  Colorado  and 
one  in  Los  Angeles,  California. 

The  first  study  was  conducted  by  Wertheimer  and  Leeper  in  1979.  In  this  study 
the  authors  designed  their  own  exposure  measurement  system  to  approximate  a 
measurement  of  long-term  exposure.  This  system  categorized  homes  according  to 
characteristics  of  the  distribution  systems  which  supplied  the  homes  with  electrical  power 
(e.g.  conductor  thickness  and  distance  of  conductor  from  the  home).  This  study  found 
that  there  was  over  a  two-fold  increase  in  the  risk  of  brain  cancer  among  children  living 
in  homes  classified  as  “high-current”  over  those  living  in  the  “low-currenf ’-classified 
homes  (see  Table  2).  Additionally,  the  risk  of  contracting  leukemia  in  high-current 
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homes  was  two  to  three  times  that  of  the  low-current  homes  (see  Table  1)  (Wertheimer 
and  Leeper,  1979:273-284). 

The  second  Denver  study  was  conducted  in  1988  by  Savitz,  et  al.  In  addition  to 
using  the  Wertheimer  and  Leeper  wiring  system,  Savitz,  et  al.,  also  recorded  short-term 
spot  measurements  of  the  electric  and  magnetic  fields  in  various  rooms  throughout  the 
homes.  As  Table  2  shows,  when  wiring  classifications  were  used  this  study  also  showed 
a  doubling  of  the  brain  cancer  risk  in  the  high-current  homes  when  compared  to  the  low- 
current  homes.  The  risk  of  leukemia,  however,  was  noticeably  lower  than  the  relative 
risks  calculated  in  the  first  (Wertheimer  and  Leeper,  1979:273-284)  Denver  study  (Table 
3)  (Savitz,  et  al.,  1988:21-38). 

The  Los  Angeles  childhood  leukemia  study  conducted  in  1991  by  London,  et  al., 
added  an  additional  measurement  to  the  second  Denver  study:  24-hour  magnetic  field 
measurements  in  the  children’s  bedrooms  (London,  et  al.,  1991:923).  This  study,  like  the 
second  Denver  study,  also  showed  a  weaker  association  between  the  wiring 
configurations  and  childhood  leukemia  than  the  first  Denver  study.  What  the  Los 
Angeles  study  did  show,  along  with  the  second  Denver  study,  is  that  there  was  no 
“consensus”  among  the  various  magnetic  field  measurements  which  would  lead  one  to 
conclude  that  there  was  an  association  with  cases  of  childhood  leukemia  (see  Table  4). 
One  implication  which  crosses  all  three  studies  is  very  important:  the  wiring 
configuration  categories  are  not  valid  indicators  of  absolute  levels  of  magnetic-field 
exposure  (Trichopoulus,  1992:V-8). 
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Table  2  RR  of  Childhood  Brain  Cancer  in  Relation  to  Wiring  Configuration,  High- 


vs.  Low-Current  Classifications  (Case-Control) 


Study 

Residence  Occupied 

Relative 

Risk 

95% 

Confidence  Interval 

Wertheimer  and 

At  time  of  birth 

2.4 

1. 1-5.1 

Leeper 1979 

At  time  of  death 

2.4 

1. 0-5.5 

Savitz,  et  ah,  1988 

At  time  of  diagnosis 

2.0 

1. 1-3.8 

2  years  before  diag. 

(Trichopoulus,  1992;V-3) 


Table  3  RR  of  Childhood  Leukemia  in  Relation  to  Wiring  Configuration,  High-  vs. 


Low-Current  Classifications  (Case-Control) 


Location 

Residence  Occupied 

i 

1 

Relative 

Risk 

95% 

Confidence  Interval 

Denver 

At  time  of  death 

3.0 

1.8-5.0 

At  time  of  birth 

'2.3 

1.3-3.9 

Denver 

At  time  of  diag. 

1.5 

2  years  before  diag. 

1.7 

0.8-3.9 

Los  Angeles 

Longest  in  etiologic 
period 

1.7 

1. 1-2.5 

(Trichopoulus,  1992:V-4) 


Table  4  RR  of  Childhood  Leukemia  in  Relation  to  Four-Category  Wiring 
Configurations  (Case-Control) 


Denver 

Los  Angeles 

Los  Angeles 

Spot 

Measurement 

Spot 

Measurement 

24-Hr 

Measurement 

Measure 

& 

Category 

RR 

95% 

Cl 

Measure  & 
Category 

RR 

95% 

Cl 

Measure  & 
Category 

RR 

95% 

Cl 

<0.65 

mG 

1.0 

<0.67  mG 

1.0 

<0.68  mG 

1.0 

0.65-1.00 

mG 

0.9 

0.3-2.7 

0.67-1.24 

mG 

1.4 

0.7- 

2.9 

0.68-1.18 

mG 

0.7 

0.4- 1.2 

1.00-2.49 

mG 

1.4 

0.5-3.4 

1.25  mG 

1.2 

0.5- 

2.8 

1.19-2.67 

mG 

0.9 

0.5-1.7 

>2.5  mG 

2.1 

0.5-7.0 

>2.68  mG 

1.5 

0.7-3.3 

(Trichopoulus,  1992:V-4) 


In  1982  and  1987  Wertheimer  and  Leeper  repeated  their  childhood  experiment  on 


adult  subjects  in  the  Denver  area  using  the  same  methodology.  They  found  that  for  all 
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cancers  combined,  the  relative  risk  for  high-current  homes  as  compared  to  low-current 
homes  was  1.4  with  a  confidence  interval  of  1.2-1. 6  (Wertheimer  and  Leeper,  1982:345- 
355;  Wertheimer  and  Leeper,  1987:43-53).  Most  of  the  other  well-known  adult  studies 
which  have  been  done  (i.e.  McDowall,  1986;  Severson,  et  al.,  1988;  Youngson,  et  al., 
1991;  Meijers,  et  al.,  1991),  have  either  been  accomplished  with  controversial  methods  or 
have  not  shown  statistically  significant  associations  between  cancers  and  wiring 
configurations. 

3.2  Occupational  Epidemiological  Research 

There  have  been  many  studies  of  occupational  exposure  to  ELF-EMFs  which  have 
sought  associations  with  various  forms  of  cancer  (i.e.  leukemia,  brain  cancer,  male  breast 
cancer,  etc.)  (Trichopoulus,  1992:V18-V22).  There  have  been  so  many,  in  fact,  that  it 
would  not  be  practical  to  list  them  all  here.  The  two  forms  of  cancer  which  were  the 
focus  of  most  of  these  studies  are  the  same  as  in  the  residential  studies,  i.e.  leukemia  and 
brain  cancer. 

In  order  to  better  view  the  overall  research  involving  risks  and  exposure, 
Trichopoulus  attempted  to  combine  the  findings  of  no  less  than  30  different 
epidemiological  studies  from  1963  to  1991  involving  cancer  and  exposure  to  EMFs. 

These  experiments  included  proportionate,  case-control,  and  prospective/retrospective 
cohort  studies  (Trichopoulus,  1992:V29-V49).  Their  method  of  aggregating  the 
experiments  consisted  of  “estimat[ing]  a  summary  proportionate  mortality  or  incidence 
ratio  for  proportionate  studies,  a  summary  odds  ratio  for  case-control  studies,  and  a 
summary  standardized  mortality  or  incidence  ratio  for  cohort  or  retrospective  cohort 
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studies”  (Trichopoulus,  1992:V-18).  The  specific  methods  used  by  Trichopoulus  to 
calculate  these  summaries  were  as  follows; 

“For  proportionate  and  cohort  studies,  the  observed  and 
expected  number  of  cases  were  obtained  first  for  each  study.  The 
summary  estimate  was  then  calculated  as  the  ratio  between  the  sum 
of  the  observed  cases  and  the  sum  of  the  expected  cases.  The 
aggregated  estimate  was  therefore  weighted  proportionally  to  the 
size  of  the  population.  The  95%  confidence  interval  for  the 
summary  estimate  was  obtained  from  a  Poisson  distribution 
(Rothman  and  Boice  1982). 

For  case-control  studies,  the  summary  odds  ratio  estimate 
was  derived  from  the  weighted  average  of  the  logarithm  of  the 
odds  ratio  of  the  individual  studies.  The  weights  were  taken  to  be 
proportional  to  the  inverse  variances  of  the  log  odds  ratio.  These 
variances  were  estimated  from  comparison  of  the  upper  and  lower 
95%  confidence  limits  on  the  assumption  that  the  difference 
between  log  upper  and  log  lower  confidence  limits  was  3.92 
standard  error  of  the  log  odds  ratio.  We  calculated  the  95% 
confidence  interval  of  the  summary  odds  ratio  by  taking  the  inverse 
sum  of  the  weight  as  variance  of  the  log  summary  odds  ratio” 

(Trichopoulus,  1992:  V- 19). 

The  results  of  this  composite  of  experiments  can  be  found  in  Table  5  and  Table  6  below. 


Table  5  Summary  Analysis  of  Occupational  Studies  on  EMFs  and  Leukemia 


Leukemia  Type 

Design 

(n  =  observed  cases) 

Summary  RR 
(95%  Cl) 

All  leukemia 

Proportional  mortality  (n  =  618) 

1.2(1. 1-1.2) 

Proportional  Incidence  (n  =  148) 

1.2(1.0-1.4) 

Case-control 

1.1  (0.9- 1.3) 

Cohort  (n  =  599) 

1.1  (1.0-1.2) 

(Trichopoulus,  1992:V-20) 


As  one  can  see  from  Table  5,  the  relative  risks  for  all  leukemia  ranged  from  1.1  to 
1.2,  however,  the  only  confidence  interval  which  could  assure  a  relative  risk  above  1.0 
was  that  for  the  proportional  mortality  studies.  These  studies  also  yielded  the  “tightest” 
confidence  interval  among  the  four  types. 
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Referring  to  the  brain  cancer  summary  table  (Table  6),  the  case-control  studies 
seem  to  show  a  fairly  significant  relative  risk  even  considering  the  confidence  interval. 
According  to  Trichopoulus,  however,  the  hypothesis  that  these  case-control  studies  were 
homogeneous  had  to  be  rejected.  The  proportional  mortality  summary  relative  risk  also 
seemed  to  be  a  bit  elevated,  but  it  was  only  based  on  one  study  and  still  included  1 .0 
within  its  95%  confidence  interval.  Finally,  the  summary  of  the  cohort  studies  did  not 
show  any  tendency  that  the  relative  risk  was  different  than  1 .0. 

What  is  interesting  to  note  is  that  for  both  leukemia  and  brain  cancer  the  widest 
confidence  intervals  are  attached  to  the  case-control  studies.  This  would  seem  to  indicate 
that  one  of  the  other  types  of  studies  would  probably  prove  to  be  of  more  statistical  value 
when  planning  new  experiments. 

Table  6  Summary  Analysis  of  Occupational  Studies  on  EMFs  and  Brain  Cancer 


Design 

(n  =  observed  cases) 

Summary  RR 

(95%  Cl) 

Proportional  mortality  (n  =  110) 

1.2 

(1.0- 1.4) 

Case-control 

1.4 

(1.2- 1.6) 

Cohort  (n  =  429) 

1.0 

(0.9-1. 1) 

(Trichopoulus,  1992:V-21) 


One  occupational  study,  however,  needs  some  specific  discussion.  This  is  the 
1994  study  by  London,  et  al.,  into  leukemia  risks  encountered  by  electrical  workers  in  the 
Los  Angeles  area  (London,  et  al.,  1994:47-60).  The  significance  of  this  study  is  that  it 
generated  and  analyzed  the  same  data  that  is  the  subject  of  this  research.  The  difference 
lies  in  the  methods  of  analysis. 

The  subjects  of  this  study  agreed  to  wear  magnetic  field  measuring/recording 
devices  which  would  measure  the  rms  magnetic  field  magnitude  every  2.5  seconds  of  one 
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work  shift.  These  values  were  then  either  averaged  or  surveyed  for  the  percentage  of  time 


that  an  individual’s  exposure  was  over  a  predefined  threshold  value.  The  resulting  values 
were  used  to  create  three  exposure  groups  (London,  et  al.,  1994:50).  The  table  of  interest 
is  Table  7. 


Table  7  OR  for  Leukemia  According  to  Estimates  of  Average  Magnetic  Field  and 
Percent  of  Workday  Above  2.5  and  25  mG  by  Occupation 


Variable  and 
Category 

Cases 

Controls 

OR  (95%  Cl) 
categorical 

OR  (95%  Cl)  per  10 
unit  increase 

Average  Magnetic  Field  mG 

<1.7 

2264 

65160 

1.0 

1.8 -8.0 

61 

1408 

1.2(1.0-1.6) 

8.1 

30 

644 

1.4(1. 0-2.0) 

All 

2355 

67212 

1.2(1.0-1.5) 

Percent  of  time>2.5  mG 

>13.0 

2264 

65160 

1.0 

13.0-32.9 

55 

1131 

1.4(1. 1-1.8) 

33.0 

36 

921 

1. 2(0.9- 1.6) 

All 

2355 

67212 

1.1  (1.0-1.2) 

Percent  of  time>25  mG 

<0.5 

2268 

65266 

1.0 

0.5  -  7.9 

57 

1302 

1.2(1.0-1.6) 

>7.9 

30 

644 

1.4(1. 0-2.0) 

All 

2355 

67212 

1.4(1. 0-2.1) 

(London,  et  al.,  1994:54) 


This  research  showed  a  trend  that  electrical  workers  were  at  a  slightly  higher  risk 
than  non  electrical  workers,  but  because  of  the  width  of  the  confidence  intervals  the 
statistical  association  was  not  very  strong. 

3.3  Laboratory  Research  Supporting  Cell  Response  Theories  of  Transiently 
Augmented  Transcription 

Significant  laboratory  research  has  also  been  conducted  which  reports  noticeable 
effects  on  various  biological  properties  as  a  result  of  externally  applied  EMFs.  First, 
certain  low-frequency,  sinusoidal  EMFs  increase  mRNA  transcription  rates  (Litovitz,  et 


al.,  1990:297).  Second,  when  such  a  sinusoidal  external  field  is  applied  and  subsequently 
remains  constant,  the  mRNA  response  rises  to  some  maximum  value  and  then  decays  to  a 
steady  state  value  which  depends  upon  the  strength  of  the  exogenous  field  (Litovitz,  et  al., 
1990:298).  Additionally,  the  magnitude  and  rise  time  of  this  maximum  response  have 
been  shown  to  depend  upon  the  strength  of  the  applied  field.  As  the  strength  of  the  field 
increases,  the  maximum  of  the  response  increases  and  the  maximum  occurs  earlier  in 
time.  This  response,  however,  does  not  increase  without  bound.  Eventually  a  saturation 
point  is  reached  and  the  maximum  response  will  no  longer  increase  with  field  strength 
(Litovitz,  et  al.,  1990:298).  Finally,  supporting  research  has  found  that  if  such  a  field  is 
prematurely  terminated  or  its  amplitude  is  modified,  the  return  to  a  steady-state  level  can 
be  greatly  affected  (Litovitz,  et  al.,  1990:305-306;  Lin,  et  al.,  1995:1).  These 
characteristics  are  modeled  in  Figure  1  and  Figure  2. 
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Figure  1  Three  Litovitz  Model  Transient  Responses 

(Litovitz,  1990:302) 

There  is  an  additional  property  which  is  not  necessarily  obvious.  This  involves 
the  fact  that  it  is  critically  important  to  know  at  what  point  in  time  the  response 
measurement  is  made  when  considering  different  field  strengths.  Looking  closely  at 
Figure  1,  as  time  progresses  the  response  to  the  strongest  field  actually  falls  below  the 
responses  of  both  the  intermediate  and  weak  fields  before  all  three  reach  a  lower  steady- 
state  value.  This  “moving”  amplitude  window  makes  timing  critical  when  taking 
measurements  of  biological  responses  to  electromagnetic  fields.  If,  for  example,  a 
researcher  were  to  record  the  affected  biological  activity  at  only  the  4-hour  mark,  he 
might  conclude,  and  incorrectly  so,  that  the  weakest  field  produces  the  largest  response. 
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If  data  were  recorded  at  the  90  minute  mark,  the  medium  strength  field  would  be 
identified,  again  incorrectly,  as  the  one  producing  the  maximum  response. 


0  2  4  6  8  10 


Time 

Unmodified  Response 
Basal  Level  for  Unmodified  Response 
”  Field  Modification  @  t=0.4 
“  Basal  Level  for  Modified  Field 
Field  Termination  @  t=0.4 
Basal  Level  for  Terminated  Field 

Figure  2  Effects  on  Transient  Response  of  Field  Variation/Termination 

Responsiveness  of  mRNA  transcription  has  been  inferred  qualitatively  from 
observed  systemic  responses.  For  example,  repetitive  pulse  train  and  repetitive  single 
pulse  EMFs  of  varying  waveforms  are  used  extensively  to  aid  healing  of  stubborn  bone 
fractures  and  in  the  treatment  of  avascular  necrosis  and  osteoporosis  (Goodman,  et  al., 
1983: 1283).  To  explore  this  directly  at  a  cellular  level,  Goodman,  et  al.,  tested  whether 
such  EMFs  had  any  effect  on  the  normal  RNA  transcription  patterns  in  the  salivary  gland 
chromosomes  of  the  dipteran  Sciara  coprophila.  Exposure  to  both  types  of 
quasirectangular,  asymmetrical,  pulsed  EMFs  resulted  in  increased  transcription  activity. 
However,  the  “effects  of  both  pulses  declined... after  60  minutes  of  continuous  exposure” 
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(Goodman,  et  al.,  1983:1283-1285)  showing  that  the  resulting  augmented  transcriptions 
were  transient. 

In  1986  Goodman  and  Henderson  expanded  this  research  to  explore  the  effects  of 
symmetrical  sine  waves  on  the  transcriptional  activity  of  Sciara  salivary  glands 
(Goodman  and  Henderson,  1986:23-29).  Sine  wave  frequencies  of  72,  222,  and  4,400  Hz 
were  applied  at  magnetic  field  strengths  of  1 1,500,  3700,  and  180  mG  respectively.  The 
72  Hz  signal  was  specifically  chosen  to  compare  with  a  single  pulsed  signal  having  a 
repetition  rate  of  72  Hz  similar  to  that  used  in  the  1983  research.  The  experiments 
revealed  that  exposure  to  the  72  Hz  sine  wave  increased  transcription  in  the  same  size 
classes  as  those  seen  following  exposure  to  72  Hz  single-pulsed  signal.  The  222  and 
4,400  Hz  frequencies  produced  the  same  general  pattern  of  induced  transcription,  but  not 
as  noticeably  as  the  72  Hz  field.  In  fact,  the  transcriptional  activity  showed  an  inverse 
correlation  with  the  field  frequencies.  One  possible  explanation  for  this  could  have  been 
the  significantly  lower  magnetic  field  strengths  of  the  222  and  4,400  Hz  frequencies, 
however  the  object  of  the  research  was  to  show  how  signal  shape,  and  not  necessarily 
strength,  would  affect  cell  transcription.  The  “results  demonstrate [d]  that  electromagnetic 
fields  within  a  specific  frequency  range,  irrespective  of  wave  shape,  can  affect 
transcriptional  activity  in  cells”  (Goodman  and  Henderson,  1986:28).  More  specifically, 
certain  symmetrical  sine  waves  as  well  as  pulsed  EMF  signals  induce  increased 
transcription  in  cells  transiently  under  a  continuous  exposure. 

The  next  logical  step  was  to  verify  that  this  augmented  transcription  could  be 
duplicated  in  human  cell  experiments.  In  1989  Goodman,  et  al.,  exposed  human 
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leukemia  HL60  cells  to  five  different  ELF  electromagnetic  signals;  three  pulsed 
asymmetric  signals  with  different  repetition  rates  (1.5,  15,  72  Hz)  and  two  continuous 
sinusoidal  waveforms  of  60  and  72  Hz  (Goodman,  et  al.,  1989b:  216-220).  The  pulsed 
signals  were  included  to  compare  with  the  earlier  Sciara  transcription  results,  the  60  Hz 
sine  wave  was  added  to  approximate  the  environment  created  by  the  power  system 
frequency,  and  the  72  Hz  sine  wave  was  included  for  comparison  with  its  corresponding 
pulsed,  asymmetric  signal.  The  result  in  each  experiment  indicated  that  the  highest 
transcription  augmentation  occurred  when  the  HL60  cells  were  exposed  to  the  sinusoidal 
signals,  and  of  those  two  signals  the  60  Hz  field  produced  the  greatest  transcription 
activity  (Goodman,  et  al.,  1989b:  218).  Additionally,  an  in-depth  analysis  of  this 
research  as  a  whole  indicated  that  within  groups  of  differing  signal  types,  the  magnitude 
of  transcript  response  seemed  to  show  a  dependence  upon  field  strength  (Goodman,  et  al., 
1989a;  Goodman,  et  al.,  1989b:2 16-220). 

The  dependence  upon  field  strength  along  with  exposure  time  was  specifically 
addressed  by  Goodman,  et  al.,  in  1992.  This  experiment,  which  exposed  human  HL-60 
cells  to  continuous  sine  waves  of  60  Hz,  found  measurable  increases  in  some  transcript 
levels  (Goodman,  et  al.,  1992: 19).  They  also  found  that  the  transcript  levels  which  were 
affected  began  responding  to  the  exogenous  fields  at  times  as  early  as  four  minutes  and 
peaked  at  twenty  minutes  when  cells  were  exposed  to  magnetic  fields  of  57  mG. 
Transcript  levels  decreased  to  near  control  levels  at  twenty  minutes  when  exposed  to 
magnetic  fields  of  5700  mG  (Goodman,  et  al.,  1992:23).  These  findings  closely  follow 
characteristics  outlined  at  the  beginning  of  this  section.  In  this  experiment  the  response 
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to  the  57  mG  field  showed  a  transient  increase  which  was  dependent  upon  the  field 
strength.  Although  not  specifically  addressed,  the  explanation  for  the  seenndngly 
unresponsive  5700  mG  field  also  supports  a  transient  response  dependent  upon  field 
strength.  The  reason  that  the  5700  mG  field  decreased  to  control  levels  after  twenty 
minutes  follows  from  the  fact  that  the  stronger  the  field,  the  earlier  the  peak  response 
occurs  and  the  quicker  the  return  to  basal  levels. 

The  most  recent  research  is  reported  in  a  1995  research  paper  by  H.  Lin,  et  al., 
currently  pending  publication  in  Bioelectrochemistrv  and  Bioenergetics.  The  aim  of  this 
research  was  to  determine  the  effects  of  both  continuous  and  limited  duration  60  Hz,  80 
mG  field  exposures  on  HL60  mRNA  transcription  activity  (c-myc)  (Lin,  H.,  et  al.,  1995: 
1-17).  The  results  showed  that  under  a  continuous  exposure  the  c-myc  transcription 
levels  peaked  after  20  minutes  and  returned  to  a  lower  steady-state  (basal)  level  after  60 
minutes.  When  the  field  was  turned  off  after  a  20  minute  exposure,  however,  it  took 
three  times  longer  to  drop  to  a  steady-state  level  (Lin,  H.,  et  al.,  1995:1).  The  significance 
of  this  research  is  that  it  was  the  first  research  which  explicitly  exhibited  the  aggregate 
properties  outlined  in  the  Litovitz  multistage  model.  This  model  will  be  covered  in 
section  3.5. 

As  noted,  these  properties  were  specifically  targeted  by  Litovitz,  et  al.,  in  1990 
when  they  proposed  a  multistage  mathematical  model  linking  magnetic  field  strength  to 
mRNA  transcription.  Development  and  application  of  this  model  are  the  subjects  of  the 
next  section.  The  adjustment  of  this  model  to  represent  the  empirical  results  of  Lin,  et  al., 
1995  is,  in  part,  the  subject  of  this  thesis. 
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3.4  Laboratory  Research  Supporting  Effects  of  Field  Coherence  on  Transiently 
Augmented  Transcription 

Although  augmented  transcription  has  proven  to  be  a  result  of  ELF-EMF 
exposure,  the  specific  exogenous  signal  characteristics  which  cause  this  to  occur  and  the 
resulting  cell  response  mechanisms  which  directly  elevate  transcription  activity  are  not 
entirely  known.  One  signal  characteristic  which  has  been  proposed  to  affect  a  cellular 
response  is  coherence.  Qualitatively,  coherence  is  inversely  related  to  the  rate  of  change 
of  amplitude  or  frequency.  A  perfectly  coherent  signal  would  have  constant  amplitude 
and  frequency.  Both  amplitude  and  frequency  coherence  have  been  studied  as  factors  in 
cellular  response. 

Byus,  et  al.,  in  1987  investigated  the  ability  of  low-energy  60-Hz  EM  fields  to 
alter  the  activity  of  ornithine  decarboxylase  (ODC)  which  is  an  enzyme  involved  in  cell 
growth.  Although  the  intent  of  this  research  was  not  to  specifically  identify  coherence  as 
the  agent  responsible  for  the  altered  enzyme  activity,  it  did  serve  as  the  springboard  for 
the  research  that  did  target  coherence.  The  results  showed  that  “a  1-hour  exposure  to  a 
60-Hz  EM  field  of  an  intensity  of  10  mV/cm  produced  a  5-fold  increase  in  ODC  activity 
in  human  lymphoma  CEM  cells  and  a  2-  to  3-fold  increase  in  mouse  myeloma  cells 
relative  to  the  unexposed  cultures.  [Further,]  depending  on  the  cell  type,  ODC  activity 
increased  during  the  1-hour  exposure  period  and  remained  elevated  for  several  hours  after 
the  field  exposure  ended”  (Byus,  et  al.,  1987:1385). 

In  1991,  Litovitz,  et  al.,  conducted  additional  experiments  using  ODC  response  as 
the  indicator.  This  time,  however,  the  intent  was  to  identify  whether  or  not  a  certain  field 
characteristic  might  be  responsible  for  increased  transcriptional  activity.  Realizing  that  a 
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cell’s  incoherent  endogenous  thermal  field  is  several  orders  of  magnitude  larger  than 
those  of  most  coherent  exogenous  fields,  Litovitz,  et  al.,  set  out  to  determine  whether  or 
not  a  cell  has  the  ability  to  discriminate  between  coherent  and  incoherent  fields.  They 
found  that  by  varying  the  frequency  between  55  and  65  Hz  with  at  least  10  seconds  at 
each  frequency  setting  (coherence  time),  using  a  100  mG  field,  and  exposing  the  cell 
cultures  for  4  hours,  they  were  able  to  enhance  ODC  activity  just  as  much  as  if  the  field 
had  been  left  at  55  or  65  Hz  for  the  full  4  hour  period.  If  a  coherence  time  of  5  seconds 
was  used,  ODC  activity  was  at  a  level  half  way  between  the  controls  and  the  10  second 
samples.  If  any  coherence  time  smaller  than  1  second  (essentially  incoherent)  was  used, 
no  activity  was  detected  (Litovitz,  et  al.,  1991:864). 

The  implications  of  these  results  were  tested  in  two  additional  papers:  Litovitz,  et 
al.,  1994a,  and  Litovitz,  et  al.,  1994b.  In  the  first  paper  Litovitz,  et  al.,  tested  whether  or 
not  the  combination  of  an  incoherent  noise  field  and  a  coherent  60-Hz  field  could  keep 
developing  chick  embryos  from  experiencing  abnormalities  which  in  earlier  experiments 
occurred  as  the  result  of  exposure  to  certain  60  Hz  fields  alone  (Litovitz,  et  al., 

1994a:  105-113).  In  the  second  paper  Litovitz,  et  al.,  explored  whether  or  not  enhanced 
ODC  activity  in  L292  cell  cultures  could  be  diminished  under  similar  conditions  as  the 
chick  embryos  (Litovitz,  et  al.,  1994b: 399-409).  In  both  instances  as  long  as  the 
superimposed  noise  was  comparable  in  magnitude  with  the  EMF  signal,  the  signal  effects 
on  the  biological  subjects  could  be  eliminated  (Litovitz,  et  al.,  1994a:  1 10;  Litovitz,  et  al., 
1994b:407). 
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3.5  Multistage  Toxicokinetic  Model 


The  definition  of  the  term  toxicokinetic  literally  means  the  movement  of  poisons 
throughout  a  body.  Although  electromagnetic  fields  are  not  technically  considered  to  be 
toxins,  the  resulting  augmentation  of  mRNA  transcription,  if  destructive,  can  be 
considered  a  toxin.  It  is  in  this  way  that  a  model  linking  EMFs  and  harmful  side  effects 
can  be  considered  toxicokinetic. 

In  1990  and  1992  Litovitz,  et  al.,  found  that  a  multistage  model  produced  an 
acceptable  approximation  of  the  properties  that  the  laboratory  data  was  revealing.  More 
specifically,  this  was  a  set  of  sequential,  first-order  chemical  reactions  modeled  by 
differential  equations.  In  order  to  posit  this  model,  however,  a  hypothesis  was  needed  to 
explain  the  multistage  mechanisms  by  which  augmented  transcription  was  achieved 
(Litovitz,  et  al.,  1990:297-312).  The  multistage  model  incorporated  with  the  chosen 
mechanisms  can  be  visualized  as  follows: 


Nucleotide 

reservoir 


Diffusion 


Positioned 
G,  C,  A,  U 


Polymerization 


Degradation 


Figure  3  Stages  of  the  Litovitz  Multistage  Model 

(Litovitz,et  al.,  1990:299) 


This  is  simplified  for  discussion: 

■■■ 

Figure  4  Variables  Representing  the  Litovitz  Multistage  Model 

(Litovitz,  et  al.,  1990:299) 

The  nucleotide  reservoir  concentration,  which  is  assumed  to  be  limitless,  is  given 
by  [A].  ki  is  the  rate  constant  representing  the  diffusion-controlled  migration  to  and  the 
positioning  and  orienting  of  the  various  nucleotides  from  the  cellular  pool  (described  in 
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chapter  2  as  A,  G,  C,  and  U).  [X]  is  the  concentration  of  positioned  and  oriented 
nucleotides  preparing  for  transcription.  k2  is  the  polymerization  rate  (oriented  nucleotides 
combining  to  form  a  transcripted  mRNA  chain)  of  the  positioned  and  oriented  G,  C,  A, 
and  U  nucleotides.  The  resulting  response,  or  mRNA  concentration,  before  degradation 
is  given  by  [Y].  And  finally,  ks  is  the  degradation  rate  of  messenger  RNA  by  cytoplasmic 
nucleases  (enzymes  that  promote  hydrolysis  of  nucleic  acids,  i.e.  chemical  process  of 
decomposition)  (Litovitz,  et  al.,  1990:299-300). 

This  multi-step  model  would,  in  reality,  consist  of  a  number  of  complex 
individual  processes,  but  the  assumption  here  is  that  the  inclusion  of  more  steps  would 
not  change  the  outcome  of  the  model  and  would  only  minimally  affect  the  shapes  of  the 
response  functions  at  the  different  EMF  strengths  (Litovitz,  et  al.,  1990:299).  The  pivotal 
hypothesis  in  the  model  is  that  switching  on  the  electromagnetic  field  at  t  =  0  produces  a 
“sudden”  increase  in  k2  and  ka.  The  increase  in  k2  and  ka  will  be  symbolized  by  k2*  and 
ka*  (Litovitz,  et  al.,  1992:241). 

If  the  assumption  is  made  that  [A],  the  nucleotide  reservoir  concentration  will 
always  remain  high  enough  to  supply  the  process,  the  first  order  sequence  of  reactions  can 
be  described  by  the  following  set  of  linear  first  order  differential  equations: 

—  =  -l-k,A-k2X  (2) 

dt  ’  ^ 


^  =  +k2X-k3y 


(Litovitz,  et  al.,  1990:300) 


(3) 
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These  values,  however  were  not  derived.  They  were  picked  by  Litovitz  to  approximate 
known  qualitative  results.  Figure  1  shows  representative  responses  for  B  =  0.2,  1.0,  and 
5.0  mG  respectively. 

There  are  two  characteristics  of  this  model  which  are  not  readily  apparent  from 
the  graphs.  One  is  the  fact  that  the  response  does  yield  a  maximum  as  a  function  of  field 
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strength.  In  other  words,  there  is  a  point  where  the  response  no  longer  increases  as  a 
result  of  field  strength  and  any  increase  in  the  field  strength-calculated  parameters  k2*  and 
ks  will  not  further  increase  the  concentration  of  mRNA  (Litovitz,  et  al.,  1992:243).  The 
peak  will,  however,  occur  earlier  as  the  field  strength  is  increased.  The  other  interesting 
behavior  is  that  a  certain  critical  time  seems  to  exist  which  will  cause  significantly  larger 
bio-effects  than  those  of  earlier  or  longer  exposure  times.  For  the  constants  used  in 
Figure  1,  the  critical  time  is  approximately  one  hour  (Litovitz,  et  al.,  1992:243). 

This  model  also  shows  that  quenching  a  field  before  the  response  has  had  a 
chance  to  settle  at  a  basal  level  has  a  significant  effect  on  the  response.  When  the  field  is 
switched  off,  the  response  not  only  takes  much  longer  to  stabilize,  it  also  tends  to 
“undershoot”  what  will  be  its  final  steady-state  value  and  then  gradually  increases  to  meet 
it  (Litovitz,  et  al.,  1990:307).  When  the  field  is  switched  off  ( t*  t ),  the  response 
function  becomes: 


y(t  >  t*)  = 

k. 


(7) 


where  the  initial  conditions  are  the  values  of  x  and  y  immediately  before  the  field  is 

terminated; 


Ak,  =  k,  -k,;  B  =  -^ 


k,A  Ak; 

kj-kj  k/ 


(l-e-^^‘  ); 


and  C  =  -B  -  ^ 

k,  k-t  J 


^3  kj 


kj-k^ 


(Litovitz,  et  al.,  1992:241-242) 
This  function  is  shown  concurrently  with  others  in  Figure  2. 
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3.6  Recent  Application  of  Multistage  Toxicokinetic  Model 

In  1994  research,  which  at  the  time  of  this  writing  is  awaiting  publication, 
Thomas,  et  al.,  reanalyzed  data  from  a  1991  case-control  study  of  childhood  leukemia  and 
electromagnetic  fields  in  Los  Angeles  County  (London,  et  al.,  1991:923-937).  They  used 
the  Litovitz  multistage  model  to  derive  an  exposure  metric  known  as  the  “kinetic  index” 
(the  integral  of  the  mRNA  response  over  the  course  of  the  measurement  period)  for  each 
child  and  then,  using  a  variety  of  statistical  analyses,  investigated  the  correlation  between 
this  index  and  the  incidents  of  childhood  leukemia  (Thomas,  et  al.,  1994: 1). 

Magnetic  field  measurements  (root  mean  squared  magnitude)  were  taken  at  fifty- 
second  intervals  in  308  residences  yielding  164  cases  and  144  controls  (Thomas,  et  al., 
1994:5).  Using  Litovitz’ s  representative  values  for  the  rate  constants,  ki,  ki,  ks,  bka,  and 
5k3,  Thomas,  et  al.,  calculated  a  response  function  for  each  separate  fifty-second  interval 
using  the  final  response  value  of  the  previous  interval  as  the  initial  response  value  for  the 
next.  The  kinetic  index  was  then  calculated  by  integrating  over  all  of  the  interval 
response  functions  of  a  particular  individual.  This  yielded  a  single  kinetic  index  metric 
per  child  (Thomas,  et  al.,  1994:8). 

The  most  notable  result  of  this  research  in  the  context  of  this  thesis  was  that 
individually  neither  the  calculated  kinetic  index  nor  the  50-second  autocorrelation  (a 
measure  in  many  regression  applications  indicating  when  the  error  terms  in  time-series 
data  are  not  independent)  trend  tests  provided  a  statistically  significant  link  to  the 
childhood  leukemia  incidents.  However,  in  a  multiple  logistic  regression  analysis  where 
both  factors  were  considered  predictors,  the  individual  trend  tests  for  the  kinetic  index 
and  the  fifty-second  autocorrelation  showed  a  statistical  link  (Thomas,  et  al.,  1994:22-23). 
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4.  Toxicokinetic  Model  Modifications  and  Application 

This  chapter  will  accomplish  three  goals.  First  it  will  outline  the  ideas  which 
were  considered  as  modifications  to  the  Litovitz  Multistage  Model.  It  will  then  chronicle 
the  logic  and  rationale  behind  the  chosen  modifications  to  the  Litovitz  Multistage  Model. 
Finally,  the  model  will  be  developed  and  used  to  create  kinetic  index  (Thomas,  et  al., 
1994: 1)  values  for  each  individual  whose  magnetic  field  measurements  were  taken  in  the 
London,  et  al.,  1994  study  (London,  et  al.,  1994:47-60). 

4. 1  Coherence-Kinetic  Modei 

Because  of  the  strength  of  association  between  the  50  second  autocorrelation  and 
the  kinetic  index  in  the  Thomas,  et  al.,  research,  the  idea  to  somehow  combine  the  two 
ideas  into  one  model  seemed  logical.  One  linking  mechanism  investigated  in  this 
research  was  posited  by  Dr.  Joseph  Bowman  of  the  National  Institute  of  Occupational 
Safety  and  Health  (Bowman,  1995). 

In  1991  Litovitz,  et  al.,  presented  a  coherence  theory  which  he  thought  could  help 
explain  augmented  cell  transcription  by  magnetic  fields.  Realizing  that  endogenous 
thermal  noise  fields  within  cells  themselves  are  much  larger  than  most  fields  caused  by 
exogenous  sources,  Litovitz,  et  al.,  tested  whether  or  not  cells  have  the  ability  to 
discriminate  between  incoherent  thermal  noise  and  weaker  coherent  (constant  phase, 
frequency,  and  amplitude)  exogenous  fields  (Litovitz,  et  al.,  1991:863). 

The  test  exposed  cultures  of  murine  L929  cells  to  60  Hz  magnetic  fields  of  10, 

100,  or  1000  mG  for  times  ranging  from  1  to  8  hours.  The  response  of  interest  was  the 
enhancement  in  activity  of  the  enzyme  ornithine  decarboxylase  (ODC)  in  the  L929  cells. 
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This  was  measured  in  terms  of  the  ratio  of  exposed  to  control  activity  (Litovitz,  et  ah, 
1991:864). 

The  link  between  augmented  mRNA  transcription  and  increased  enzyme  activity 
is  direct.  Enzymes  are  constructed  from  proteins.  The  rate  of  synthesis  of  a  protein  is 
proportional  to  the  amount  of  its  corresponding  mRNA  that  is  present,  any  increase  in  the 
amount  of  mRNA  will  subsequently  result  in  an  enhanced  rate  of  protein  production  and, 
therefore,  enzyme  production. 

The  results  of  the  Litovitz,  et  al.,  experiment  showed  that  “application  of  [constant 
amplitude]  fields  for  four  hours  but  with  [frequency]  coherence  times  of  10  or  50  seconds 

[produced]  enhancements  in  ODC  activities . In  contrast,  for  coherence  times  of  0. 1  or 

1 .0  seconds  [essentially  incoherent]  no  enhancement  of  ODC  activity  was  observed” 
(Litovitz,  et  al.,  1991:864).  Further  results  showed  that  “for  the  cell  to  respond  to  an  ELF 
signal  it  is  necessary  for  the  exogenous  field  to  maintain  coherence  for  a  minimum  time 
interval  greater  than  about  [5]  seconds,  with  full  response  requiring  an  interval  greater 
than  about  10  seconds”  (Litovitz,  et  al.,  1991:864). 

In  order  to  further  explore  these  observations,  in  1994  Litovitz,  et  al.,  conducted 
experiments  to  see  whether  or  not  superimposing  temporally  incoherent  (both  frequency 
and  amplitude)  magnetic  fields  onto  coherent  60  Hz  fields  suppressed  the  ODC  activity  of 
exposed  L929  cells  (Litovitz,  et  al.,  1994a:399-409).  “It  [was]  concluded  that  the 
superposition  of  incoherent  magnetic  fields  can  block  the  enhancement  of  ODC  activity 
by  a  coherent  magnetic  field  if  the  strength,  or  rms  amplitude,  of  the  incoherent  field  is 
equal  to  or  greater  than  that  of  the  coherent  field”  (Litovitz,  et  al.,  1994a:399).  By  curve- 
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fitting  data  from  this  experiment,  Litovitz  developed  the  following  equation  which 
describes  ODC  activity  as  a  function  of  the  noise-to-signal  ratio  (a  measure  of  amplitude 
coherence): 


[ODC]  =  l  + 


1.06 


1  +  76 


N 

vS  y 


(8) 


where  N  =  noise  (rms)  and  S  =  signal  (rms). 

(Litovitz,  et  al.,  1994a:405) 

Figure  5  shows  the  relationship  of  the  ODC  activity  ratio  with  the  noise-to-signal  ratio. 
The  smaller  the  noise-to-signal  ratio,  the  greater  the  amplitude  coherence  and  the  greater 
the  resulting  ODC  activity  will  be.  Unfortunately,  this  equation  does  not  include  a  factor 
for  frequency  coherence,  noted  in  Litovitz,  et  al.,  1991  as  an  additional  factor  affecting 
ODC  activity. 
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ODC  Activity  Function 


Figure  5  ODC  Activity  Ratio  vs.  Noise-to-Signal  Ratio 

(Litovitz,  et  al.,  1994a:406) 

As  a  result  of  the  multistage  model  and  the  experiments  relating  coherence  to 
enzyme  activity,  Dr.  Bowman  proposed  combining  the  two  models  as  the  subject  of  this 
thesis.  Originally  this  was  to  be  accomplished  by  modifying  the  multistage  rate  constants 
using  a  factor  which  is  controlled  not  only  by  the  magnetic  field  strength,  i.e. 
kj*  =  kj  +  B  •  (Skj)  hr  '  and  kj*  =  kj  +  B  •  (bkj)  hr'',  but  also  by  the  coherence,  i.e. 

kj*  =  kj  +  B  (5k2)  C  hr  '  and  =  kj  +  B  (5k3)  C  hr  ',  with  coherence  (C)  being 
unitless  (Bowman,  1995). 

The  data  used  in  the  London,  et  al.,  research  in  1994  was  taken  by  magnetic  field 
meters  which  were  worn  by  individuals  representing  various  occupations  for  one  duty 
day.  The  meters  recorded  magnetic  field  rms  values  (B)  every  2.5  seconds  (London,  et 
al.,  1994:49-50).  In  order  to  calculate  a  coherence  factor  from  equation  8,  a  means  of 
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computing  noise  and  signal  had  to  be  established.  Bowman  proposed  estimating  the 
signal  term  by  the  moving  average  rms  field  value  over  the  previous  ten  second  period 
(four  2.5-second  intervals).  He  did  so  because  he  interpreted  the  results  of  Litovitz,  et  al., 
1991  to  show  that  the  minimum  time  for  which  coherence  seemed  to  show  noticeable 


effects  on  ODC  activity  to  be  5  to  10  seconds  (Litovitz,  et  al.,  1991:864).  The  noise  term 
would  then  be  the  difference  between  the  moving  average  and  the  actual  measured  rms 
field  value  during  the  specific  2.5-second  time  increment  of  interest.  Equations  9,  10,  and 
1 1  show  specifically  how  these  values  were  calculated. 


noise,  =  11.  -SIGNAL,]' 


(Litovitz,  et  al.,  1994a:405) 

Each  COHERENCEi  term  is  represented  by  C  in  the  linear  rate  constant  equations  that 
include  coherence. 


There  were,  however,  a  few  of  problems  with  this  idea  for  a  kinetic-coherence 
model.  The  first  problem  surfaced  in  Bowman’s  interpretation  of  the  Litovitz,  et  al.,  1991 
results.  Litovitz  did  not  observe  increased  ODC  activity  after  10  seconds  of  exposure  to 
an  amplitude  and  frequency  coherent  field.  What  his  results  showed  was  that  if  an 
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amplitude  coherent  field  varied  its  frequency  in  10-second  intervals  or  longer,  e.g. 
oscillating  between  55  Hz  for  the  first  10  seconds  and  65  for  the  next  10-second 
increment,  for  a  period  of  four  hours,  maximum  ODC  activity  would  be  observed  after 
this  four-hour  period. 

The  second  problem  lay  in  the  1994  research  by  Litovitz,  et  al.  This  experiment 
was  partially  discussed  earlier  in  section  4.1  and  resulted  in  equation  8.  Even  though  this 
equation  is  a  fairly  good  predictor  of  the  ODC  activity  ratio,  the  noise  that  was 
superimposed  on  the  60  Hz  field  not  only  was  amplitude  incoherent,  it  also  varied 
frequency  between  30  and  90  Hz.  Granted,  noise  frequency  fluctuations  are  more 
difficult  to  analyze,  but  since  Litovitz,  et  al.,  had  demonstrated  in  1991  that  this  frequency 
variation  could  not  be  ignored,  equation  8  should  have  reflected  some  sort  of  frequency 
dependence. 

Even  if  Litovitz,  et  al.,  1994  had  measured  and  shown  a  frequency  dependence, 
the  data  used  in  this  thesis  (London,  et  al.,  1994)  contained  no  frequency  information 
other  than  the  fact  that  the  B -field  readings  were  within  the  40-400  Hz  frequency  range. 
With  this  in  mind,  frequency  fluctuations  could  not  be  used  to  develop  a  modified 
multistage  response  model. 

The  final  result,  not  really  a  problem,  which  was  the  decisive  factor  in  discarding 
the  idea  to  combine  the  kinetic  index  and  coherence  combination  was  revealed  in  a 
modeling  comparison.  When  a  sample  data  set  for  a  single  individual  was  used  in  both  a 
model  which  included  Bowman’s  suggested  coherence  term  and  a  model  which  was  the 
result  of  this  research  and  did  not  include  coherence,  there  were  essentially  no  major 
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differences  in  model  output.  Figure  6  and  Figure  7  below  show  that  although  the  mRNA 
responses  are  not  completely  identical,  they  do  indicate  that  the  coherence  term,  when 
used  in  conjunction  with  smaller  time  constants,  derived  in  this  research,  plays  nearly  no 
role  when  compared  to  a  similar  response  model  using  only  the  smaller  time  constants. 

mRNA  Response  with  Coherence 


Time  (seconds) 

mRNA  Measurement  x  10 
■  ■  Actual  B  Field  Measurements 

Figure  6  Modeled  mRNA  Response  with  Coherence 
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mRNA  Response  Without  Coherence 


Actual  B  Field  Measurements 


Figure  7  Modeled  mRNA  Response  without  Coherence 

Since  this  specific  combination  of  the  coherence  and  kinetic  models  did  not 
appear  to  be  proper,  from  this  point  on  coherence  was  not  considered  to  be  part  of  the 
model. 

4.2  Modified  Muitistage  Modei 

During  the  research  of  this  thesis  Lin,  et  al.,  released  results  which  seemed  to 
directly  support  the  characteristics  demonstrated  by  the  Litovitz  multistage  model.  Their 
laboratory  data  was  recorded  under  conditions  of  constant  phase  and  frequency  (60  Hz.) 
with  the  only  varying  factor  being  the  amplitude  (Lin,  et  al.,  1995).  This  research 
reported  complete  response  times  which  were  on  the  order  of  one  to  two  hours. 

Since  there  was  no  laboratory  data  which  specifically  supported  the  numerical 
values  of  the  multistage  rate  constants’  chosen  by  Litovitz,  et  al.,  1991  and  1992,  the 


38 


focus  of  this  thesis  became  the  derivation  of  new  rate  constants  for  the  Litovitz  multistage 
model  combined  with  an  application  of  the  resulting  model  to  investigate  leukemia 
incidents  in  certain  Los  Angeles  County  electrical  occupations. 

Using  equation  4,  normalized  so  that  the  steady-state  ambient  field  response  was 
1,  and  the  data  points  which  were  supplied  by  Lin,  et  al.,  1995,  we  attempted  to  obtain 
values  for  ki,  k.2,  k2*,  ks,  and  ks*.  We  found  that  by  limiting  the  algorithm  to  only  four 
degrees  of  freedom  we  could  not  obtain  an  adequate  fit  of  a  new  function  to  the  data 
points.  This  is  clearly  shown  in  Figure  8. 


Response  Model  with  Constant  kl 


Time  (hours) 

Model  Graph  for  Continuous  Field 
■  ■  Model  Graph  for  Terminated  Field 
Data  Plot  for  Continuous  Field 
Data  Plot  for  Terminated  Field 

Figure  8  Response  Model  Using  Lin,  et  al.,  Data  and  Constant  ki 

An  additional  indicator  suggested  that  this  was  an  unrealistic  fit  of  the  data  as 
well.  The  value  for  the  rate  constant  ki,  which  at  this  point  had  no  corresponding  6ki 
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value,  was  an  extremely  large  negative  number. 


v98 


,  -5.196  10’'  .1 

k,  = - sec  . 

'  60 


The  other  rate 


constants  were  on  the  order  of  10”^  to  10^  and  were  positive.  Even  given  the 
speculative  nature  of  the  underlying  physical  model,  this  value  of  ki  was  completely 
unrealistic. 

At  this  point,  an  additional  question  was  added  to  the  focus  of  the  research.  Could 
ki,  in  addition  to  k2  and  ka,  also  be  accelerated  by  the  presence  of  a  magnetic  field?  This 
was  a  very  logical  question  since  both  k2  and  ks  both  had  corresponding  6k  values  By 
adding  this  new  parameter,  we  were  able  to  obtain  an  acceptable  fit  of  the  resulting 
function  to  the  data  points.  The  calculations  of  the  new  parameters  are  found  in 
Appendices  G  and  H.  The  plots  of  the  data  and  resulting  function  are  found  in  Figure  9. 
The  parameters  are: 


,  2.357 -10-" 

k,  = - sec 

’  60 


,  1.965  10^ 

k,  = - sec 

"  60 


kj  = 


-1.383  IQ-^ 
60 


■  sec 


and  the  rate  constants  for  any  interval,  i 


kj  *  =  kj  +Bi  -bk, , 


^2,  = 


6k,  = 
bk^  = 

bkj  = 


1.214  10~^ 
60 

7.152  10^ 
60 

7.361  10^ 
60 


sec''mG’’ 


sec'mG'* 


sec''mG'* 


=  1  to  n,  would  be: 

kj  +  Bj  -bkj,  and  kj  *  =  k3  +  Bj  •  bkj 
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Data  Points  &  Resulting  Function 


Time  (hours) 

Continuous  Function 
■  “  Function  Terminated  at  20  min. 

Data  Points  Under  Cont.  Field 
Data  Points  of  Field  Term.  @  20  min. 


Figure  9  Response  Model  Using  Lin,  et  al.,  Data  and  Variable  ki 

It  is  interesting  to  note  how  much  improvement  was  obtained  by  adding  an 
additional  degree  of  freedom  to  the  model.  One  would  intuitively  expect  improvement  by 
adding  a  new  parameter,  but  this  improvement  should  be  relative  to  the  number  of 
degrees  of  freedom.  In  this  case  the  model  has  shown  an  unusually  strong  affinity  to  the 
additional  parameter. 

Based  upon  the  original  Litovitz  model  ki  “represents  the  diffusion-controlled 
migration  to,  and  the  positioning  and  orienting  of  the  various  nucleotides  from  the 
cellular  pool”  (Litovitz,  et  al.,  1990:299)  along  the  DNA  strands  in  the  nucleus  of  the  cell. 
One  key  assumption  in  the  Litovitz  model  is  that  this  part  of  the  multistage  process  is 
passive  and  is  unaffected  by  the  presence  of  a  magnetic  field.  However,  as  noted  earlier. 
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this  was  a  speculative  assumption  and  is  likely  to  be  an  oversimplification  of  the 
underlying  biological  processes.  It  is  plausible,  then,  that  ki  may  actually  be  affected  by 
the  presence  of  a  magnetic  field  just  as  k2  and  ks  are  theorized  to  be. 

On  the  other  hand,  this  is  still  highly  speculative.  We  know  only  that  if  a 
multistage  model  such  as  this  is  appropriate,  then  it  must  include  variation  of  ki  with  the 
magnetic  field.  We  also  know  that  this  multistage  model  exhibits  all  of  the  observed 
qualitative  phenomena  and  fits  the  known  quantitative  data.  The  underlying  mechanisms, 
however,  are  still  unclear. 

4.3  Application  Procedure 

The  next  step  was  to  use  the  new  rate  constants  in  the  Litovitz  multistage  model  to 
calculate  the  magnetic  field  responses,  or  kinetic  indexes  (Thomas,  et  al.,  1994:1),  for 
each  2.5  second  interval  for  each  individual  in  the  London,  et  al.,  1994  study.  This  was 
accomplished  using  a  Matlab®  program  (see  Appendix  A  for  source  code)  incorporating  a 
finite  difference  algorithm  (see  Appendix  B).  The  output  from  this  program  will  be 
analyzed  in  the  next  chapter. 
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5.  Data  Analysis 

This  chapter  should  enable  the  reader  to  understand  how  the  output  of  the 
computer  simulated  model  was  calculated  from  the  magnetic  field  data  taken  from  the 
workers  in  Los  Angeles.  Additionally,  one  should  comprehend  the  rationale  and 
statistical  bases  behind  the  chosen  statistical  groupings.  Next,  odds  ratio  calculations, 
which  were  the  ultimate  result  of  the  statistical  analysis,  should  allow  the  reader  to 
understand  the  conclusions  which  will  be  outlined  in  the  next  chapter.  Finally,  one 
should  be  able  to  distinguish  the  differences  between  using  the  average  magnetic  field 
and  the  kinetic  index  as  indicators  of  exposure  on  identical  data  sets. 

5. 1  Kinetic  Index  Calculations 

Although  there  has  been  some  question  during  this  research  regarding  the  proper 
rate  constants  which  should  be  used  when  applying  the  Litovitz  multistage  model,  the 
model  itself,  as  we  have  shown,  is  flexible  enough  to  accurately  represent  the  targeted 
biological  properties  with  only  slight  modifications.  This  and  the  fact  that  the  rate  of 
change  of  the  response  function  is  almost  negligible  when  considering  only  2.5  second 
intervals  led  to  a  greatly  simplified  method  of  calculating  kinetic  indexes. 

The  following  linear  first  order  differential  equations,  (2)  and  (3),  were  first 
presented  in  chapter  three  (Section  3.5).  Instead  of  using  the  solution  to  these  equations 
to  calculate  the  incremental  responses  from  the  data  set,  the  responses  were  approximated 
at  each  interval  by  a  set  of  finite  difference  equations  (12)  and  (13). 
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(2) 


—  =  +k,A-k,x 
dt  ‘  ' 

dy 

dt  ' 

(Litovitz,  et  al.,  1990:300) 

XX, =xxi  +dt-(k,(b;)-xx.  •k2(b.))  (12) 

and  yyi.,  =yyi+dt-(k2(bi)-xxi-k3(bi)-yyi)  (13) 

where  ki(bi),  k2(bi),  and  k3(bi)  are  the  “starred”  k  values  from  section  3.5,  each  bi  value  is 
the  magnetic  field  reading  found  in  the  i'**  recorded  data  interval,  and  dt  =  2.5  seconds  is 
the  length  of  each  interval.  Once  each  yy;  value  was  calculated,  all  were  totaled  for  each 
individual  and  multiplied  by  2.5  seconds  in  order  to  get  an  approximation  for  the  time 
integral  over  the  duty  day  (kinetic  index).  These  individual  indexes  were  then  averaged 
over  job  codes  representing  the  different  occupations  of  the  monitored  individuals.  The 
occupational  kinetic  indexes  along  with  means,  sample  sizes,  and  standard  deviations  are 
found  in  Table  8. 
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Table  8  Mean  Kinetic  Index  and  Associated  Statistics  by  Occupation 


Occupation 

Mean  K.I. 
(“dose”-niinutes)  ' 

Sample  Size 

Std.  Dev. 

“Exposure 

Group 

Nonelectrical  Worker 

410.74 

95 

144.50 

Electrical  Engineer 

439.41 

14 

179.17 

Low 

Phone  Line  Worker  and 
Splicer 

450.84 

21 

139.96 

Electrician  and 
Apprentice 

491.62 

28 

168.44 

Motion  Picture 
Projectionist 

497.79 

13 

241.19 

Medium 

Electrical  Engineering 
Technician 

507.97 

11 

128.30 

TV  and  Radio 
Repairman 

568.28 

19 

173.73 

Electric  Power  Wire 
and  Cable  Worker 

585.43 

87 

148.40 

High 

Welder  and  Flame 
Cutter 

656.48 

22 

132.81 

Power  Station  Operator 

686.76 

34 

135.49 

Total 

523.22 

344 

152.95 

In  order  to  show  a  comparison  with  the  kinetic  index  we  also  calculated  the 
average  magnetic  field  exposure  for  each  occupation.  Even  though  this  was  the  identical 
exposure  indicator  and  data  used  by  London,  et  al.,  1994,  it  was  not  possible  to  use  their 
figures  directly  because  they  task-weighted  the  averages.  Since  the  report  did  not  indicate 
how  this  weighting  was  accomplished,  the  averages  had  tabe  recalculated  directly  from 
the  data.  The  occupational  means,  sample  sizes,  and  standard  deviations  are  found  in 
Table  9. 
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Table  9  Mean  Average  Magnetic  Field  Exposure  and  Associated  Statistics  by 
Occupation 


Occupation 

Mean  Average 
Magnetic  Field  (mG) 

Sample  Size 

Std.  Dev. 

“Exposure” 

Group 

Nonelectrical  Worker 

1.5081 

95 

1.0498 

Low 

Electrical  Engineer 

1.8478 

14 

1.8819 

Phone  Line  Worker 
and  Splicer 

1.9046 

21 

1.1355 

Electrician  and 
Apprentice 

2.91 

28 

1.9976 

Medium 

Electrical  Engineering 
Technician 

3.0778 

11 

2.3403 

TV  and  Radio 
Repairman 

3.9836 

19 

2.6887 

Motion  Picture 
Projectionist 

7.918 

13 

7.0589 

High 

Power  Station 
Operator 

13.942 

34 

25.952 

Welder  and  Flame 
Cutter 

17.059 

22 

24.376 

Electric  Power  Wire 
and  Cable  Worker 

32.047 

87 

70.157 

Total 

12.036 

344 

152.95 

5.2  Statistical  Groupings,  Tests,  and  Odds  Ratios 

In  1994  London,  et  al.,  “performed  a  registry-based  case-control  study  among  men 
aged  20-64  years  with  known  occupations  who  were  diagnosed  with  cancer  in  Los 
Angeles  County  between  1972  and  1990”  (London,  et  al.,  4994:47).  The  case-control 
values  from  this  study  are  found  in  Table  10.  The  raw  data  from  the  monitored  subjects 
with  which  the  case-control  data  were  to  be  compared  are  the  same  used  in  this  research. 
London,  et  al.,  however,  used  the  average  magnetic  field  as  the  “treatments”  as  opposed 
to  the  modified  kinetic  index  derived  in  this  research. 
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Table  10  Cases  and  Controls  by  Category  of  Electrical  Occupation 
Among  Men  Aged  20-64  and  Diagnosed  with  Cancer  in  L.A.  County, 
1972-1990 


Occupational  Category 

Leukemias 

(Cases) 

Controls 

Nonelectrical  Worker 

2234 

64547 

Low 

Electrical  Engineer 

30 

613 

Phone  Line  Worker  and 
Splicer 

4 

31 

Electrician  and  Apprentice 

28 

728 

Medium 

Motion  Picture  Projectionist 

1 

22 

Electrical  Engineering 
Technician 

24 

i 

521 

1 

TV  and  Radio  Repairman 

4 

106 

High 

Electric  Power  Wire  and 
Cable  Worker 

2 

50 

Welder  and  Flame  Cutter 

27 

579 

Power  Station  Operator 

1 

15 

Total 

2355 

67212 

(London,  et  al.,  1994:54) 


Table  11  Mean  KJ.  and  Associated  Statistics  by  Occupational  Grouping 


Group 

Mean  KI 
(“Dose”-min.) 

95%  Cl 
(“Dose”-min.) 

Sample  Size 

Group  Std. 
Dev. 

Low 

420.30 

394.72  -  445.89 

130 

147.46 

Medium 

496.62 

446.84  -  546.41 

52 

178.82 

614.34 

590.68  -  638.00 

162 

152.49 

As  with  London,  et  al.,  we  decided  to  use  three  groupings  of  occupational 
categories  with  the  difference  that  they  were  grouped  by  average  kinetic  index  rather  than 
average  magnetic  field  exposure.  These  groups  were  denoted  as  low,  medium,  and  high 
“exposure.”  An  attempt  was  made  to  ensure  that  each  group  would  include  at  least  thirty 
cases  of  leukemia.  This  was  done  so  that  there  would  be  enough  cases  of  leukemia  to 
make  statistically  valid  comparisons  between  the  chosen  groups.  These  groupings  are 
found  in  both  Table  8  and  Table  10. 
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In  order  to  assure  that  the  three  groupings  were  statistically  different  and  had 
similar  variances,  two  tests  were  accomplished.  The  first  test  was  Bartlett’s  test  of  equal 
variances  (see  Appendix  C).  In  this  test  the  null  hypothesis  posited  that  the  variances  of 
all  three  groups  were  equal.  The  alternate  hypothesis  stated  that  not  all  variances  were 
equal.  With  our  kinetic  index  groupings,  the  p-value  for  the  test  was  0.2245  indicating 
that  only  if  the  significance  level  were  greater  than  this  could  we  reject  the  null 
hypothesis.  We  could,  therefore,  not  reject  the  null  hypothesis  that  all  three  variances 
were  equal. 

The  second  test  was  Tukey’s  pairwise  comparison  of  means  which  would  indicate 
whether  or  not  the  three  groups  had  means  which  were  significantly  different  from  one 
another.  Upon  execution  of  the  test  we  found  that  all  three  means  were  significantly 
different  from  one  another  at  a  rejection  level  of  0.050  (see  Appendix  C). 

The  final  steps  consisted  of  calculating  the  case/control  ratios,  the  subsequent 
odds  ratios  and  the  confidence  intervals  for  the  odds  ratios  for  all  three  groups  (see 
Appendix  E  for  calculations).  Odds  ratios  and  confidence  intervals  are  shown  in  Table 
12. 


Table  12  K.I.  Odds  Ratios  and  C.Ls  by  Occupational  Grouping 


Group 

Odds  Ratios 

Confidence  Intervals 

Low  Dose 

1 

N/A 

Medium  Dose 

1.199 

0.9081  -  1.5821 

High  Dose 

1.303 

0.9227-  1.8403 

Alternatively,  we  also  grouped  occupations  by  average  magnetic  field  exposure. 
This  was  the  exposure  metric  used  by  London,  et  al.,  1994,  and  the  one  which  has  not 
given  consistent  results.  Again  we  chose  three  occupational  groupings  (low,  medium. 
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and  high)  with  an  emphasis  on  including  at  least  thirty  leukemia  cases  in  each  group. 
These  groupings  are  found  in  Table  9.  Even  though  this  was  the  best  choice  of  several 
tested  groupings,  there  were  no  low-medium-high  groupings  of  occupations  where 
Bartlett’s  test  of  equal  variances  was  not  rejected.  Essentially,  this  meant  that  at  least  two 
groups  of  any  chosen  three  based  upon  average  occupational  magnetic  field  exposure 
would  not  have  variances  which  were  statistically  indistinguishable  (see  Appendix  D). 

There  was  also  no  combination  which  led  to  the  finding  of  significantly  different 
means  between  all  three  chosen  groups  (Appendix  D).  This  was  the  most  disconcerting 
finding  in  that  if  groups  cannot  be  statistically  separated,  then  any  statistical  calculations 
based  upon  them  become  meaningless.  Statistics  on  the  chosen  groupings  based  upon 
average  occupational  magnetic  exposure  can  be  found  in  Table  13. 


Table  13  Mean  Magnetic  Field  Exposure  and  Associated  Statistics  by 
Occupational  Grouping 


Group 

Mean 
Magnetic 
Exp.  (mG) 

95%  Cl 
(mG) 

Sample  Size 

Group  Std.  Dev. 

Low 

1.6088 

1.4043  -  1.8132 

130 

1.1781 

Medium 

3.2935 

2.6840  -  3.9030 

58 

2.3181 

23.977 

15.248  -  32.705 

143 

55.188 

As  with  the  kinetic  index  groupings,  the  final  step  Was  to  calculate  case/control 
ratios,  odds  ratios,  and  the  confidence  intervals  for  these  odds  ratios  (see  Appendix  F). 
Odds  ratios  and  confidence  intervals  are  shown  in  Table  14. 


Table  14  Avg.  Mag.  Field  Exposure  Odds  Ratios  and  C.I.s  by 
Occupational  Grouping 


Group 

Odds  Ratios 

Confidence  Intervals 

Low  Dose 

1 

N/A 

Medium  Dose 

1.23 

0.9433  -  1.6047 

High  Dose 

1.338 

0.9322  -  1.9202 
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One  will  immediately  notice  that  the  odds  ratios  of  the  medium  and  high  groups 
of  the  average  magnetic  field  measurement  groupings  are  slightly  higher  than  those  of  the 
kinetic  index  groupings  (Table  12).  The  groupings  for  the  average  magnetic  field 
measure,  however,  are  questionable  due  to  the  fact  that  statistically  separate  occupational 
groupings  could  not  be  constructed  using  average  magnetic  field  as  the  measure.  In  this 
case  the  inability  to  create  statistically  separate  groupings  may  lead  one  to  conclude,  and 
not  incorrectly  so,  that  odds  ratio  calculations  based  on  average  magnetic  field  rather  than 
the  kinetic  index  are  statistically  invalid. 

An  interesting  comparison  exists  between  the  two  measurement  methods’  odds 
ratio  confidence  intervals.  For  both  the  medium  and  high  groupings  the  separate 
measurement  method  confidence  intervals  have  essentially  the  same  widths.  Here  the 
explanation  is  not  so  much  in  the  validity  of  the  groupings  as  it  is  in  the  quality  and 
numerical  characteristics  of  the  case-control  data.  Even  though  we  compared  two 
completely  different  measurement  methods  for  EMF  exposure,  the  same  case-control 
numerical  data  was  used  for  correlation  purposes.  Per  grouping  we  did  not  significantly 
change  the  numerical  values  enough,  between  measurement  methods,  to  adversely  affect 
the  widths  of  the  confidence  intervals. 

Finally,  we  note  the  correlation  between  the  kinetic  index  and  the  arithmetic  mean 
of  the  magnetic  field  exposure  for  all  of  the  monitored  subjects.  One  will  recall  that 
Thomas,  et  al.,  1994  used  Litovitz’s  rate  constant  values  and  did  not  allow  ki  to  vary  in 
the  analysis  of  its  Los  Angeles  childhood  leukemia  data.  In  the  Thomas  research  the 
correlation  was  almost  perfectly  negative  (p=-0.97)  using  the  log  values  (Thomas,  et  al.. 
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1994:25).  The  correlation  using  the  London,  et  al.,  1994  data  set,  allowing  ki  to  vary,  and 
using  the  rate  constants  derived  in  this  research  showed  a  nearly  completely  opposite 
trend  (p=0.64),  also  using  log  values.  This  positive  correlation  makes  physical  sense, 
since  we  would  expect  that,  if  there  is  a  biological  response,  generally  greater  magnetic 
fields  should  give  generally  elevated  responses. 
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6.  Conclusions 


This  final  chapter  should  tie  all  aspects  of  this  research  together.  By  now  the 
reader  should  have  a  good  grasp  on  the  rationale  behind  modifying  Litovitz’s  multistage 
model  as  well  as  an  understanding  of  how  the  changes  were  made.  The  conclusions  and 
recommendations  will  be  discussed  in  the  areas  of  model  development  and  model 
application  separately. 

6. 1  The  Improved  Toxicokinetic  Model 

The  most  beneficial  result  of  this  research  was  the  improvement  of  the  Litovitz 
multistage  model.  Recent  experimental  data  supplied  by  Lin,  et  al.,  1995  allowed  the 
Litovitz  multistage  model  parameters  to  be  calculated  directly  from  biological  data.  Until 
this  time,  no  calculation  of  parameters  had  been  accomplished  directly  from  biological 
data. 

Of  similar  importance  was  the  discovery  that  the  biological  data  could  be 
accurately  modeled  only  by  allowing  all  three  of  the  time  constants  to  vary  with  the 
imposed  field.  This  required  adding  an  additional  parameter,  5ki.  The  assumption  by 
Litovitz,  et  al.,  that  ki,  described  as  diffusion-controlled  positioning  of  the  nucleotides 
(Litovitz,  et  al.,  1990:300),  does  not  change  under  magnetic  field  exposure  was 
inconsistent  with  the  Lin,  et  al.,  data.  This  could  mean  one  of  two  things.  Either  the 
biological  processes  which  Litovitz  chose  to  explain  the  separate  stages  of  the  multistage 
model  are  not  sufficient,  which  Litovitz,  et  al.,  have  implied  may  be  the  case  (Litovitz,  et 
al.,  1990:299),  or  diffusion  controlled  positioning  of  the  nucleotides  is  indeed  affected  in 
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some  manner  by  exposure  to  magnetic  fields.  More  research  will  be  needed  to  determine 
the  proper  conclusion. 

Another  suggestion  for  additional  research  would  be  another  experiment  to 
establish  more  data  points  for  the  transient  behavior  of  mRNA  concentration  along  the 
same  lines  as  Lin,  et  al.,  1995.  Although  there  were  sufficient  data  points  to  calculate  the 
modified  model  parameters  in  this  research,  more  data  points  would  increase  the 
confidence  in  the  accuracy  of  the  fit  of  the  model.  In  particular,  it  would  be  useful  to 
quantify  the  response  in  the  first  few  minutes  to  help  insure  accurate  modeling  of 
changing  fields. 

6.2  Model  Application 

One  interesting  difference  between  the  kinetic  index  and  the  average  magnetic 
field  measures,  is  that  slightly  different  occupational  groupings  were  established  by  each. 
Taking  a  look  at  Table  8  and  Table  9  one  will  find  that  the  low  groupings  are  identical  in 
their  occupational  compositions.  The  medium  exposure  groups  both  contain  electricians 
and  electrical  engineering  technicians,  but  the  third  occupation  is  different.  The  average 
magnetic  field  exposure  table  (Table  9)  contains  the  motion  picture  projectionist  while 
the  kinetic  index  exposure  table  (Table  8)  moves  this  occupation  to  the  high  exposure 
group  and  replaces  it  with  TV  and  radio  repairmen  from  the  average  magnetic  field’s  high 
exposure  group. 

Additionally,  the  occupations  which  are  part  of  both  high  exposure  groups,  are  not 
in  the  same  order  of  ascending  values.  The  average  magnetic  exposure  shows  power 
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station  operators  with  the  highest  exposure,  while  the  kinetic  index  yields  electric  power 
wire  and  cable  workers  as  its  highest  exposure  occupation. 

The  most  noticeable  result  of  the  application  of  the  modified  Litovitz  multistage 
model  is  the  fact  that  the  point  estimates  for  the  odds  ratios  using  the  kinetic  index  as  the 
measure  of  exposure  were  not  quite  as  high  as  those  calculated  using  average  magnetic 
field  exposure.  One  would  expect  that  they  would  be  different  due  to  the  different 
methods  of  measuring  exposure.  The  inability  to  create  statistically  separate  groupings 
using  average  magnetic  field  as  the  measurement  is  another  possibility.  Another 
explanation,  however,  is  that  the  average  magnetic  field  metric  may  inflate  the  true  odds 
ratios  by  encouraging  inappropriate  groupings. 

The  main  reason  that  the  application  of  the  model  developed  in  this  research  did 
not  yield  any  more  conclusive  results  (based  on  Cl  calculations)  than  did  London,  et  al., 
in  1994  was  due  to  the  case-control  type  of  epidemiological  study  conducted  by  London, 
et  al.  All  of  the  subjects  in  the  study  were  in  a  Los  Angeles-based  cancer  registry.  The 
cases  were  defined  as  those  subjects  who  had  leukemia.  The  controls  were  identified  as 
those  subjects  with  other  forms  of  cancer  (London,  et  al.,  1994:48). 

The  first  problem  is  that  data  begins  with  the  bias  that  all  participants  had  some 
form  of  cancer.  Additionally,  it  is  more  than  likely  that  these  records  were  somehow 
historically  incomplete  thereby  adding  an  additional  element  of  bias.  Finally,  incidence 
rates,  by  the  nature  of  a  case-control  study,  are  not  available  since  no  population-at-risk 
(total  number  of  people  who  can  develop  disease  or  population  at  the  beginning  of  a 
cohort  study)  has  been  defined.  This  effectively  eliminates  calculation  of  risk  directly  in 
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terms  of  the  most  desirable  statistic -relative  risk  (the  best  epidemiological  risk  measure) 
and  odds  ratios  become  the  best  estimates  of  relative  risk  by  default. 

Most  epidemiologists  and  statisticians  would  agree  that  either  a  prospective  or 
retrospective  cohort  study  would  have  narrowed  the  confidence  intervals  for  both  types  of 
exposure  measurements.  A  cohort  study  would  also  allow  the  calculation  of  relative  risk 
rather  than  using  the  odds  ratio  as  an  estimate.  Either  way  the  statistical  comparison  of 
the  two  measurement  methods  would  carry  more  weight.  Although  long  term  and 
expensive,  a  cohort  study  is  the  best  way  to  improve  upon  this  research. 

6.3  Closing  Comments 

The  model  developed  in  this  research  is  definitely  more  appropriate  for  modeling 
cellular  biological  responses  of  exposure  to  low  frequency  electromagnetic  fields.  This  is 
simply  because  data  existed  where  it  had  not  allowing  a  biologically  valid  model  to  be 
constructed. 

It  may  initially  seem  that  the  modified  kinetic  index  seems  to  be  no  better  at 
measuring  human  system  response  to  EMFs  than  the  average  magnetic  field.  This  is  due 
to  nearly  identical  CIs  and  lower  point  estimates  for  the  odds  ratios.  Actually,  we  believe 
it  to'be  a  more  realistic  measure  due  to  its  biological  basis.  With  this  in  mind,  the 
modified  kinetic  index  becomes  a  more  conservative  measure  of  exposure  (yielding  lower 
odds)  when  compared  to  the  average  electromagnetic  field,  and  supports  the  position  that 
ELF-EMF  exposure  does  not  create  as  great  of  a  biological  effect  as  is  currently 
hypothesized. 

It  will  take  much  more  research  to  establish  whether  or  not  the  modified  kinetic 
index  is  the  best  measure  of  ELF-EMF  exposure,  but  one  cannot  ignore  the  validity  of  its 
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biological  basis.  One  also  has  to  consider  the  fact  that  no  exposure  metric  to  date  has 
proven  conclusive  in  predicting  incidents  of  cancer.  The  fact  that  the  modified  kinetic 
index  is  a  much  more  conservative  approach  to  true  exposure  may  help  prove  its  validity 
indirectly  as  well. 
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Appendix  A.  Matlab  Source  Code  for  K.L  Calculations 

%  A  MATLAB  M-file  to  process  the  EMDEX 
%  magnetic  field  exposure  data  from  NIOSH. 

% 

%  currentfy ,  1  )=IDENTinC ATION  NUMBER 
%  current(y,2)=mG  MAGNETIC  FIELD  LEVEL 
%  current(y,3)=JOB  CODE  OF  THE  INDIVIDUAL 
%  currency ,4)=REGION  (i.e.  1=LA,  2=WA,  3=NZ) 

%  currency, 5)=k  1(b) 

%  current(y,6)=k2(b) 

%  current(y,7)=k3(b) 

%  currently, 8)=xx  concentration 
%  currently ,9)=yy  concentration 
% 

clear; 

format  long  e; 
load  mag.dat; 

%  The  original  data  file  has  been  modified  to  only  include  the  original 
%  field  values  known  as  ID,  M,  BOC,  and  REGION. 

%  The  REGION  field  was  modified  as  follows  due  to  MATLAB 's  inability  to 
%  read  ASCII  data: 

%  LA  =  1 

%  WA  =  2 

%  NZ  =  3 


limit=max(size(mag))+l ; 

%  Needed  to  add  one  line  (above)  to  the  data  file  to  compare  last  id  #  to. 


mag(limit,:)=[0  0  0  0]; 

%  For  the  above  purpose,  I  pick  an  arbitrary  array  value  which  will  not 
%  match  the  last  individual's  ID  field. 


scale=l/60; 

%  This  will  change  the  output  scale  to  dose-minutes. 


x=l;binit=4; 

%  binit  is  the  ambient  magnetic  field  level  used  in  Dr.  Reba  Goodman's 
%  experiment  (4  mG). 

kl=2.3574227589E-04*scale;delkl=1.2144280925E-03*scale; 
k2=-1.3834909426E-03*scale;delk2=7.3608201871E-04*scaie; 
k3=1.9646532847E-04*scale;  delk3=7.1523014986E-04*scale; 

%  This  sets  all  of  the  rate  constants. 
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klO=kl+binit*delkl;  k20=k2+binit*delk2;  k30=k3+binit*delk3; 

%  This  sets  the  initial  rate  constants  based  on  the  ambient  magnetic 
%  field  referenced  above. 

xx0=kl0/k20; 

yy0=kl0/k30; 

%  These  are  the  initial  values  of  the  X  and  Y  concentrations  in  the 
%  multistage  model  as  calculated  from  the  ambient  mag.  field  noted  above. 

cnt=0; 


while  x<=limit-l; 

%  Pointer  is  at  the  first  record  in  an  individual’s  file 


y=i; 

cnt=cnt+l; 

kinetic(cnt,l)=mag(x,l);  kinetic(cnt,2)=mag(x,3);  kinetic(cnt,3)=mag(x,4); 

current(y,:)=mag(x,:); 

currency, 5)=kl+current(y,2)*delkl; 

currency, 6)=k2+current(y,2)*delk2; 

current(y  ,7)=k3  +current(y  ,2)  *  delk3 ; 

currency, 8)=xx0+2.5*(kl0-xx0*k20); 

currency, 9)=yy  0+2. 5*(k20*xx0-k30*yy0); 

x=x+l; 

y=y+i; 


%  Pointer  is  at  the  second  record  in  an  individual's  file. 


beg=x; 

while  mag(x,  1  )==mag(x- 1,1); 
x=x+l; 

y=y+i; 

end; 

%x 

%y 

%beg 

current(2: 

y,l;4)=mag(beg:x,l:4); 

current(2;y,5)=kl+current(2:y,2)*delkl; 

current(2:y,6)=k2+current(2:y,2)*delk2; 

current(2;y,7)=k3+current(2:y,2)*delk3; 


for  i=2:y; 

current(i,8)=current(i- 1 ,8)+2.5*(current(i- 1 ,5)-current(i- 1 ,8)*current(i- 1 ,6)); 
current(i,9)=current(i- 1 ,9)+2.5*(current(i-  l,6)*current(i- 1 ,8)-current(i- 1 ,7)*current(i- 
1,9)); 
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end; 

kinetic  (cnt,4)=(2 . 5/60)  *  sum(current( :  ,9)) ; 

%  save  c:\school\thesis\current.dat  current  /ascii; 

clear  current; 
end; 

save  kinetic.dat  kinetic  /ascii; 
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Appendix  B.  Finite  Difference  Calculations 


Kinetic  model 

The  differential  equation  will  be  approximated  by  coupled  difference  equations  with  time 
varying  coefficients.  The  base  time  unit  will  be  one  second,  and  the  output  will  be  stored  at 
dt  second  intervals.  The  imposed  field  will  be  interpolated  between  samples. 

Pick  up  data  EMDEX  ::=READPRN(magdata)  DEX -EMDEX"'^^  length(DEX)  =  5.664- 10^ 
data  every  dt. data  seconds  dt^(^  :=2.5  -dt^j^^g  lengthtDEX) 

Set  up  rate  constants  scale  =— 

60 

k  iQ  =  2.3574227589- lO'^^  scale  k2o  =-  1.3834909426- 10'^  scale  k3o  =  1.9646532847- lO  '^-scale 

5k  j  =  1.2144280925- 10"^-scale  8k2  =7.3608201817- 10'“*  scale  8k3  :=  7.1523014986- 10'“*-scale 

k|(b)  =k2Q+b-8kj  k2(b)  :  =  k2o+ b-8k2  k3(b)  :  =  k3Q+ b-8k3 

Set  up  integration  and  storage  parameters 

to  =0  =leiigth(DEX)-2.5  Opjg  :=length(DEX)  i  :=0..np2g-l  dt  :  =  2.5  tt.  =i-dt 


b^DEX. 

^xx.  \  /  k  j(b)  -  xx.-k2(b) 

'  +  dt- 

\Pi]  \k2(b)-xx.-k3(b)-yyj 


mRNA  Response  Without  Coherence 


Time  (seconds) 

—  mRNA  Measurement  x  10 
■  ■  Actual  B  Field  Measurements 
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Appendix  C.  Kinetic  Index  Metric  Statistics 

STATISTIX  4.1®  LA_DATA,  1 1/01/95,  13:20 

ONE-WAY  AOV  FOR  KI  BY  GROUP 

SOURCE  DF _ SS _ MS _ F _ P 

BETWEEN  2  2.759E+06  1.379E-t-06  57.50  0.0000 

WITHIN  341  8.180E-1-06  23987.6 

TOTAL  343  1.094E-t-07 

« 

CHI-SO  DF  P 

BARTLETT'S  TEST 

EQUAL  VARIANCES  2.99  2  0.2245 

COCHRAN'S  Q  0.4154 

LARGEST  VAR  /  SMALLEST  VAR  1 .4705 

COMPONENT  OF  VARIANCE  FOR  BETWEEN  GROUPS  1 2863.9 

EFFECTIVE  CELL  SIZE  105.4 


GROUP 

MEAN 

SAMPLE 

SIZE 

GROUP 

STD  DEV 

1 

420.30 

130 

147.46 

2 

496.62 

52 

178.82 

3 

614.34 

162 

152.49 

TOTAL 

523.22 

344 

154.88 

CASES  INCLUDED 

344 

MISSING  CASES 

TUKEY  (HSD)  PAIRWISE  COMPARISONS  OF  MEANS  OF  KI  BY  GROUP 

HOMOGENEOUS 

GROUP  MEAN _ GROUPS 

3  614.34  I 

2  496.62  I 

1  420.30  I 

ALL  3  MEANS  ARE  SIGNIFICANTLY  DIFFERENT  FROM  ONE  ANOTHER. 

CRITICAL  Q  VALUE  3.3 14  REJECTION  LEVEL  0.050 

STANDARD  ERRORS  AND  CRITICAL  VALUES  OF  DIFFERENCES 

VARY  BETWEEN  COMPARISONS  BECAUSE  OF  UNEQUAL  SAMPLE  SIZES. 
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Appendix  D.  Average  Magnetic  Field  Metric  Statistics 

STATISTIX  4.1®  LAAVGMAG,  1 1/01/95,  13:21 

ONE-WAY  AOV  FOR  AVG_B  BY  GROUP 


SOURCE  DF 

BETWEEN  2 

WITHIN  341 

TOTAL  343 


SS _ MS 

40810.7  20405.3 

4.726E+05  1385.84 

5.134E-H05 


F _ P 

14.72  0.0000 


CHI-SO  DF  P 

BARTLETT'S  TEST 

EQUAL  VARIANCES  1080.17  2  0.0000 


COCHRAN'S  Q  0.9978 

LARGEST  VAR  /  SMALLEST  VAR  2194.6 


COMPONENT  OF  VARIANCE  FOR  BETWEEN  GROUPS  177.463 
EFFECTIVE  CELL  SIZE  1 07.2 


GROUP 

MEAN 

SAMPLE 

SIZE 

GROUP 

STD  DEV 

1 

1.6088 

130 

1.1781 

2 

3.2935 

58 

2.3181 

3 

23.977 

156 

55.188 

TOTAL 

12.036 

344 

37.227 

CASES  INCLUDED 

344 

MISSING  CASES 

TUKEY  (HSD)  PAIRWISE  COMPARISONS  OF  MEANS  OF  KI  BY  GROUP 

HOMOGENEOUS 

GROUP  MEAN _ GROUPS 

3  23.977  I 

2  3.2935  I 

1  1.6088  I 


THERE  ARE  2  GROUPS  IN  WHICH  THE  MEANS  ARE 
NOT  SIGNMCANTLY  DIFFERENT  FROM  ONE  ANOTHER. 

CRITICAL  Q  VALUE  3.3 14  REJECTION  LEVEL  0.050 
STANDARD  ERRORS  AND  CRITICAL  VALUES  OF  DIFFERENCES 
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Appendix  E.  Odds  Ratios  for  the  Kinetic  Index  Metric 


Kinetic  Index  Groupings 


Low  Exposures 

Electrical  Engineer  (30/613) 

Nonelectrical  worker  (2234/64547) 

Phone  line  worker  and  splicer  (4/3 1)  [2268/65191] 


Medium  Exposures 

Electrician  and  apprentice  (28/728) 

Motion  picture  projectionist  (1/22) 

EE  Technician  (24/521)  [53/1271] 

Hi2h  Exposures 

TV  and  radio  repairman  (4/106) 

Electric  power  wire  and  cable  worker  (2/50) 
Welder  and  flamer  cutter  (27/579) 

Power  station  operator  (1/15)  [34/750] 
CaseControl  Ratio  r  /-« 

CaseControl_Ratio 


CaseControlRatio  jyj[£D 
CaseControlRatio  j^q 


®^MED 


1.199 


CaseControlRatio 

OR  TTF  • — -  OR  Tj-r  =  1 .303 

CaseControlRatio  j^q 


CaseControlRatio 


2268 

65191 


CaseControlRatio 


53 

1271 


CaseControl  Ratio 


HI  - 


34 

750 


D 

D 

E 

A 

B 

N. 

E 

C 

D 

N. 

M, 

M, 

T 

Xsq(A-Mi,M2,Ni,N2.T) 


Nj  M  j\^ 
T  / 

N  1-N2-M  J-M2 
T^.(T-  1) 


chi mfp,  : = X  sq( 53 , 232 1 , 66462 ,1324, 67459 , 68783 ) 
chi  = X  sq( 34, 2302 , 6594 1 , 784 , 67459 ,68243 ) 


CIi^,^gr(Z,chi_square,OR)  :=OR^ 
^^lower(^-^^’‘^*“MED’®^MED)  =0.9081 
CIlower(l-96,chiHi,ORHl)  =0-9227 


^Iupper(2’‘'J“-S‘l«are,OR):=OR' 

upper( ^  06, chi  j^,  OR  =  1 .8403 
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Appendix  F.  Odds  Ratios  for  Average  Magnetic  Field  Measurements 


Average  Magnetic  Field  Groupings  „  ^  .  2268 

CaseControl  Ratio  T  n  - 

65191 


Low  Exposures 

Electrical  Engineer  (30/613) 

Nonelectrical  worker  (2234/64547) 

Phone  line  worker  and  splicer  (4/31)  [2268/65191] 

Medium  Exposures 

Electrician  and  apprentice  (28/728) 

TV  and  radio  repairman  (4/106) 

EE  Technician  {2AIS2\)  [58/1355] 


High  Exposures 

Motion  picture  projectionist  (1/22) 

Electric  power  wire  and  cable  worker  (2/50) 
Welder  and  flamer  cutter  (27/579) 

Power  station  operator  (1/15)  [31/666] 

CaseControl  Ratio  t  o 

oRlo  - -  oRlo  =  1 

CaseControI_RatioLo 


®Rmed 


CaseControlRatio 
CaseControl_Ratio  £q 


^RmED"R23 


OR  HI 


CaseControl  Ratio  hi 
CaseControlRatio  i^q 


OR  HI  ~  1.338 


CaseControlRatio  jyjgQ 


58 

1355 


CaseControl  Ratio 


HI  - 


31 

666 


D 


D 


E 

A 

B 

N. 

E 

C 

D 

N. 

M, 

M, 

T 

Xsq(A,MHM2,N 


N  I'N  2  M  j  M  2 
T^(T-  1) 


chi  : = X  sq( 58, 2326 , 66546 ,1413, 67459 , 68872 ) 

chi  HI  :=  X  sq( 5 1 , 2299, 65857 , 697 , 67459 ,68156) 


Cl  ^ chi_square,  OR)  :=  OR '  Vchi.squarc 

^^lower(^-^^’‘^^MED’^RMED)  =0-9433 
Cl lower(l-96,chi  HI, ORm)  =0.9322 


CI^ppgr(Z,chi_square,OR)  :=OR^  Vd“_square, 

CIupper(^-^^’‘=*“MED’C>RMED)  "  ^6047 
C^upper(C96,chiHi,ORHi)  =  1  9202 
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Appendix  G.  Mathcad  Parameter  Fitting  Routine 


Experimental  fit  procedure  for  the  kinetic  model. 

The  basic  model  equation  for  the  [Y]  concentration  level  is: 


f(t,XQ,yo,ki,k2,k3 


■  +  - 


^2-^3 


■exp(-k2-t)  + 


-k3-  k2-k3 


■exp/-k3t 


Note  that  there  is  no  scaling  included. 

The  corresponding  equation  for  the  [X]  concentration  level  is: 


kj  k2XQ-kj 

gT.XQ.k  J,k2)  :  =  —  + - ■exp('-k2  t 


For  the  curve  fitting  we  will  use  a  more  complicated  conditional  that  takes  both  response  curves 
and  puts  them  in  successive  time  intervals.  We  also  force  consistency  in  the  initial  conditions  this 
way.  The  scaling  to  set  ambient  steady-state  to  1  is  included. 

'  /  P4  P4  \  P6 

f  t, — , — >Pi»P2’P3  ' —  if  t <260 
\  P5  P6  /  P4 

otherwise 


fit(^t,p  i,P2.P3.P4>P5-P6) 


xO-g 

/  P4  1 

20, — ,p  2,p2 

\  P5  / 

y0<— f| 

'  P4  P4 

20,  — , — ,p  j,p  2»P  3 

(  P5  P6 

) 

f(t-  280,x0,>'0,p4,p  5,P6)- 

P6^ 

P4/ 

The  curve  fit  routine  requires  a  vector  function  that  returns  the  fijnction  value  at  time  t  plus  the  partial 
derivatives  of  the  function  with  respect  to  all  the  parameters.  Here  we'll  use  a  numerical  approximation 
for  the  derivatives. 
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D(t,p)  =  pl^Po 
p2^Pi 
P3^P2 

p4^P3 

P3^P. 

P6^P5 

fit(t,pl,p2,p3,p4,p5,p6) 

— fit(t,pl,p2,p3,p4,p5,p6) 
dpi 

—  fil(t,pl,p2,p3,p4,p5,p6) 
dp2 

—  fit(t,pl,p2,p3,p4,p5,p6) 
dp3 

- — fit(t»pl»p2,p3,p4,p5,p6) 
dp4 

— fit(t,pl,p2.p3,p4,p5,p6) 
dp5 

^ — fit(t,pl,p2.p3,p4.p5.p6) 

1  L<lp6 

Now  we  furnish  the  data  from  Lin,  et  al.,  1995  in  order  to  start  the  curve  fit. 

tdat  =(10  20  60  120  280  320  380  440  680)^ 


Xjgj  =(1.25  1.35  1.07  1.05  1.37  1.35  1.23  1.03  0.95)'^ 


The  first  four  are  the  unswitched  curve,  the  last  five  are  the  decay  part  of  the  switched 
curve. 


jl  =0..3  j2  =4..8 

Now  we  furnish  initial  guesses  for  the  parameter  values,  p  =  j—  —  —  —  —  i 

\10  15  20  200  500  350/ 


And  now  the  actual  fit. 

p  =genfit(t^t,Xjj^t,Pijji,D) 

Genfit  is  a  least  squares  curve 
fitter  using  a  gradient  search 
method  to  find  the  parameters 
that  minimize  mean  squared 
errors. 


9.7389989672-10 
5.7503070511-10 
5.7414877317-10" 
5.0934546457-10 
1.5608371301-10 
3.0573859279-  lO" 


Figure  9  shows  the  graph 
resulting  from  these  parameters. 
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Appendix  H.  Simultaneous  Equation  Solver  for  Parameters  in  Appendix  G 

In  Appendix  G  parameters  were  defined  in  a  way  which  was  easier  for  Mathcad  to  use  its  genfit 
function,  ki,  ki*,  k2,  k2*,  ks,  and  ks*  are  p4,  pi,  ps,  p2,  p6,  ps  respectively  in  Appendix  G.  In  order 
to  get  the  full  picture,  however,  kl,  k2,  k3,  6ki,  8k2,  and  fiks  need  to  be  separated  out  of  the  starred 
variables  and  defined  in  terms  of  a  4  mG  ambient  magnetic  field  which  was  used  in  Lin,  et  al.,  1995. 
Below,  a=ambient,  or  the  unstarred  variables  above,  and  e=accelerated,  or  the  starred  variables  above 

kjg  =5.0934546457- 10'^  k2a  =  1.560837130M0'^  k33  =3.0573859279-10"^ 

kjg  =9.7389989672-10"^  k2e  =5.7503070511-10"^  k3g  =5.7414877317-10"^ 

63  =4  Ambient  magnetic  field  magnitude,  Lin,  et  al.,  1995. 

B  g  =—  B  g  =80  Experimentally  applied  exogenous  magnetic  field  magnitude,  Lin,  et  al.,  1995. 

Initial  guesses  for  parameter  values  which 
can  be  inserted  directly  into  equations  4  and  5. 

Given 


kio=0  k2o  =0  k3g  =0 

8k  I  =0  8k2  =0  8k3  =0 


*'3a*‘^'3o-^^a'5^3 


•^le^kio-Bg-Skj 
‘^2e=*^2o  +  Bg-8k2 
*^30*^30+^®  e'^^3 


Find^k  jQ,8k 


2.3574227589- 10~‘* 
1.2144280925-10’'^ 
-1.3834909426-10”^ 
7.3608201817-10”“* 
1.9646532847- 10”“* 
7.1523014986-10”'* 


Final  Parameter  Values 
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It  is  a  well  known  fact  that  high  frequency,  or  radio  frequency,  radiation  can  be  directly  harmful  to  biological 
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;  radiation  to  cause  adverse  biological  effects.  Two  specific  measure,  average  magnetic  field  and  kinetic  index, 

have  been  used  in  studies  to  determine  if  they  are  useful  predictors  of  adverse  effects,  specifically  leukemia. 
Unfortunately,  both  have  shown  only  marginal  results. 

I  This  thesis  analyzed  and  improved  the  kinetic  index  measure  (Thomas,  et  al.,  1994)  by  utilizing  laboratory 

biological  data  to  increase  the  validity  of  the  parameters  of  an  existing  model.  Then,  using  this  “improved” 
model,  raw  magnetic  field  exposure  data  from  a  previous  sturfy  was  analyzed  and  compared  with  the  measurement 
method  used  in  that  sturfy  (average  magnetic  field). 

The  results  of  the  model  modification  showed  that  accurate  modeling  of  the  existing  laboratory  data  could 
only  be  accortqilished  by  including  and  additional  parameter  to  the  existing  Litovitz  multistage  model.  The 
results  of  the  model  application,  however,  were  not  particularly  conclusive.  The  case-control  study  used  for  the 
!  model  application  was  less  than  ideal,  requiring  that  more  rigorous  epidemiological  studies  be  conducted  in  order 
to  accurately  test  the  improved  model. 
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